DYNAMIC STABILITY OF EQUILIBRIUM CAPILLARY DROPS 
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Abstract. We investigate a model for contact angle motion of quasi-static capillary drops resting on a 
horizontal plane. We prove global in time existence and long time behavior (convergence to equilibrium) 
in a class of star-shaped initial data for which we show that topological changes of drops can be ruled 
out for all times. Our result applies to any drop which is initially star-shaped with respect to a a small 
ball inside the drop, given that the volume of the drop is sufficiently large. For the analysis, we combine 
geometric arguments based on the moving-plane type method with energy dissipation methods based 
on the formal gradient flow structure of the problem. 
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1. Introduction 

We consider the function u{x,t) : x [0,oo) — > [0,oo) solving the following free boundary problem 

-Au{x,t) = \{t) in {u > 0}, 

(P-V) 

V = F{\Du\) on a{w>0}, 
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Figure 1. Description of the problem 



where the Lagrange multipher X{t) is associated to the volume constraint 

J u{-,t)dx = V. (1.1) 

(see Figure 1.) 

V = Vx,t given above denotes the (outward) normal velocity of the free boundary d{u > 0} at 
{x,t). The prescribed normal velocity _F : M — >■ M is a continuous and strictly increasing function with 
F{1) — 0, where 1 denotes the equilibrium contact angle of the profile u with the surface. See Section 2 
for a rigorous formulation of weak solutions for ( P-V I . 



In two dimensions, (P-V) is a simplified model to describe a liquid droplet resting on a flat surface. 



The model is quasi-static, the speed of the contact line is much slower than the capillary relaxation time: 
see below for more discussion on the derivation of the model (also see e.g. [2], (l6]). In this context, 
u denotes the height of the drop. We call fit '■— {u{-,t) > 0} the wetted set and the free boundary 
Ft := c?{u(-,t) > 0} the contact line between the liquid and the fiat surface. The function 

F : [0, +oo) ^ M 

represents the dependence of the normal velocity of the contact line Fj on the contact angle. The 
deviation of the contact angle |Vm| on Ft (it) from the equilibrium value 1 is responsible for the motion 
of the contact line. The relation between the velocity of the contact line and the contact angle |Vit| is 
not a settled issue, and various velocity laws have been proposed and studied (e.g. [s], [l6], [Ts]) in the 
fiuids literature. This motivates us to study our problem (P-V) with the most general possible F, at 
least for the study of geometric properties of solutions. 

In section 4 and 5, we will focus on F{\Du\) — \Du\'^ — 1 for simplicity of presentation, since the 
associated energy structure is simplest for this choice of F. The approach we present nevertheless will 
apply to the general F, as we will discuss in section 4. 



Our aim is to address the long time behavior of the solution of (P-V I. Note that the equilibrium 



solution v{x) for (P-V) solves 



-Av{x) = A in {v > 0}, 
\Dv\ = 1 on d{v > 0}. 



(EQ) 



Assuming that v is C^'", the classical result of Serrin yields that v is radial. We would like to 
show the dynamic stability of the above result, in the context of our model. More precisely, we would 
like to show that the solutions of ( |P-V[ ) uniformly converges to the round drop given by (EQ|. Such 
a stability result is, to the best of the authors' knowledge, new in any model describing evolution of 
volume-preserving drops, (see [l] and [t] for stability results on stationary problems). 
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In general, the behavior of evolving drops is highly difficult to observe due to the diverse possibility of 
topological changes the drop may go through, and due to the generic non-unique nature of the evolution 
(see [2] for examples). On the other hand, it is reasonable to expect that the drops stay simply connected 
if its initial shape is star-shaped with respect to some ball inside of the drop. This is what we confirm in 



terms of p-reflection (see Definition 3.10) which measures star-shapedness of a drop with respect to the 



reflection comparison. Numerical simulations suggest that initially convex drops may develop an wedge 
while it contracts (see e.g. [15| ), and thus our definition, which allows wedges on the free boundary, 
seems to be appropriate to address the global-time behavior of the drop. 



Theorem 1.1 (Theorem 5.1). Suppose uq satisfies p-reflection and J uq = V, with < p < jq^^^^ 
Then the following holds: 



(a) there exists an "energy" solution (see the definition in section 5) u(x,t) of (P-Vl with initial 



data Uq which stays star-shaped for all times and is Holder continuous in space and time variable. 



(b) Any energy solution of (P-V) with initial data uq converges uniformly modulo translation to the 



unique equilibrium solution v of (EQ) with volume V. Moreover the set {u{-,t) > 0} converges, 



modulo translation, to a ball in Hausdorff topology. 



Remark 1.2. It should be pointed out that the result is not a perturbative one. In fact for any initial base 
of the drop {uq > 0} which is star-shaped with respect to some ball inside it (e.g. a star or a triangle), 
our theorem implies given that the drop has sufficiently large volume. Indeed the proof is based on 
geometric, moving-plane method-type arguments (see section 3.6). As a consequence we obtain explicit 
estimates on the size of the parameter p. 

Remark 1.3. It may be possible to obtain the rate of the convergence by further investigation of the 
formal energy dissipation inequality carried out in the beginning of section 4. Unfortunately in our 
setting u is not regular enough for the calculation to go through. 

Remark 1.4. Corresponding results were proved for convex solutions of the volume preserving, anisotropic 
mean curvature flow by Belletini et. al. [4], but their approach strongly depends on the level set formu- 
lation of the problem as well as the convexity preserving and regularizing feature of the mean curvature 
flow. 



It is an open question whether there exist other geometric properties besides p-reflection that are 



preserved throughout the evolution (P-V I. In fact, finding such a geometric property is one of the main 



novel features in our result. Let us point out that, in particular, it is unknown whether the convexity of 



the drop is preserved in the system ( P-V ) 



1.1. The Model. The energy of a static droplet which occupies a subset of resting on the plane 

{xN+i = 0} is given by 

J{E)^ P{Er\{xN+i>0})+ adn^ + i pgxN+idx. (1.2) 

"'En{xjv + i=0} JE 

Here P{F) is the perimeter of a subset F of defined at least for bounded domains with Lipschitz 
boundaries. Stationary droplets correspond to minimizers of the energy J when the volume \E\ — V is 
fixed. The first term is the energy due to surface tension of the free surface of the droplet. The second 
term is the energy due to the adhesion between the droplet and the surface it rests on, a is called the 
relative adhesion coefficient and in general it can be spatially dependent. In this paper however we will 
assume a is constant. The third is a gravitational potential energy, p is the mass per unit volume of 
water and g the acceleration due to gravity. We will neglect the effects of gravity and assume that 

one expects this to be a reasonable approximation when the total volume is small. 
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We will make the further assumption that our droplets are the region under the graph of a function 

u:R^ ^ [0,+oo), 

E = {{x',XN+i) e X [0,+oo) : xn+i < u{x')}. 
Now, thinking of the droplet in terms of the function u, its height profile, the energy simplifies to become, 



J{u)= ^Jl + \Du\^dx + a\{u>0}\. 



(1.3) 



M>0 



The Euler-Lagrange equation corresponding to fixed volume energy minimizers involves the mean cur- 



vature of the free surface. Finally, to simplify our analysis, we linearize (1.3) and taking a = —1/2 for 
concreteness we get 



Jiu) = 



[ \Du\^dx + \{u> 0}\. 



(1.4) 



Formally our problem (P-V) can be written as the gradient flow associated with the energy J^{u) (see the 
heuristic argument in section 4). Based on this structure a minimizing movement scheme in the spirit 
of 24 was carried out for the general energy J{E) in [2j to derive a weak solution in the continuum 
limit, where the solution is given as an evolution of entire drop (sets of finite perimeters in IR^+^). In 
our paper we show that, in the graph setting with the energy J^{u), for a large class of initial data (see 



Remark 1.2) we are able to give a point- wise description on the movement of the contact line d{u > 0} 



using the notion of viscosity solutions. 



1.2. Main challenges and strategy. The major difficulty of this problem is the lack of strong com- 
pactness (either in the associated energy given above or in the problem itself) which hinders most 
approximation arguments without geometric restrictions. In particular the standard comparison princi- 
ple fails due to A(t), or, more fundamentally, the volume preserving nature of the problem. We emphasize 
that one cannot simplify the problem by replacing A fixed in the problem, especially when one is inter- 
ested in the long-time behavior of the solutions: in fact, one can show that, if we fix A then most drops 
either vanish in finite time or grow to infinity. In comparison to other volume-preserving problems such 
as the volume-preserving mean curvature flow or the Hele-Shaw flow with surface tension, the evolution 



of our problem (P-V I is driven by the contact angle, not by higher-order regularizing terms such as the 
mean curvature. This aspect of the problem challenges, for example, the diffusive-interface approach 
which are present in many other problems. 



To get around aforemented difficulties and establish the existence of solutions for ( P-V ) , we utilize 



the energy structure of the problem. Formally, solutions of the problem ( P-V ) arc gradient flows of an 
energy associated to the capillary drop problem. This formulation of the problem is best understood 
from the perspective of the wetted set. The energy of a capillary drop of volume V resting above a 
domain C is taken to be 

J{n) := [ \Du\'' + \n\ (1.5) 
Jn 

where, for now, u is the height profile above the wetted set and is defined by 



:= arg mm{u e H^{il.) : / u = V}. 



(1.6) 



Then one models the motion of the capillary droplet as a gradient fiow of the functional in an 
appropriate space of subsets of . 

A rigorous justification of this approach is carried out using a regularized discrete gradient fiow scheme 
in the space of Cacciopoli sets (sets of finite perimeter) in [Tt] (also see [2] for more general approach). 
Since the analysis has been performed in the general setting which allows pinching and merging of droplet 
components, the resulting continuum solution is rather weak and is unstable under the variation of initial 
data. In this paper we show that much stronger results can be obtained when the initial data satisfies 
the p-refiection (see Theorem |1.1| . Our analysis will be based on the aforementioned energy structure 
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of the problem with geometric arguments (reflection maximum principle) as well as a modified viscosity 
solution theory where one views X{t) as a prescribed parameter. 

Let us explain the geometric argument in more detail. For a priori given positive and continuous 
A(t), let us consider the problem 



(P-A) 



-Au{x,t) = X{t) in {u > 0}, 

ut^ F{\Du\)\Du\ on ^{-u > 0}. 
For a given function X{t) as above and under some assumption on the niceness of the initial data fig 



and boundedness and ellipticity assumptions on F, a comparison principle holds for (P-A) and, by a 
standard application of Perron's method, there exist global-in-time viscosity solutions. 



We show, under an additional assumption on the star-shapedness of the initial data (p-reflection) , 



the existence of such a function X{t) defined for all times such that a viscosity solutions u of (P-A I 
satisfies 

u{x, t) dx — j u{x, 0) dx =: V 

for all i > 0. The proof is carried out by showing that, when fio is strongly star-shaped, then the "energy 
solution" , constructed with the discrete time scheme associated with gradient flow structure mentioned 



above, stays star-shaped and coincides with the viscosity solution of (P-A) with the corresponding A 



We point out that it is not a priori clear whether the discrete-time energy solutions preserve the star- 
shaped condition, so we incorporate the restriction directly into the approximate scheme (See section 
4) to obtain the continuum limit. The introduction of geometric restriction to the gradient flow scheme 
seems to be new and of independent interest. 



1.3. Outline of the paper. In section 2 we recall known results about the solutions of [EQ). It turns 
out that viscosity solutions of {EQ) with connected and Lipschitz positive phase are radial and minimize 
the associated energy. 



In section 3 we investigate the geometric properties of the viscosity solutions of ( P-A ) . In particular 
we show that p-reflcction property is preserved over the time (see Corollary 3.28), based on the reflection 
maximum principle. We point out that the reflection maximum principle holds only when the solution 



of (P-A) is stable under perturbations, which is the case when the solutions are star-shaped (see section 
3.5 as well as the Appendix). 



In section 4 we discuss the energy structure of our original problem ( |P- V ) . Motivated by the for 



mal gradient flow structure, we construct a solution of (P-V) as a continuum limit of a discrete-time 
"minimizing movement" scheme following the approach of 11 , 24 and [2]. By putting the p-reflection 



property as a constraint for the minimizing movements, we obtain uniform convergence of solutions in 
the continuum limit. An important result proved here is Proposition 4.7, which states that the limiting 



"energy" solution is indeed a viscosity solution of (P-A) with volume-preserving property. 



In section 5 we make use of the energy structure of the problem to show that any energy solution 
obtained in section 4 uniformly converges to the radial solution as t — >■ oo, modulo translation (The- 
orem 5.1) We point out that our result does not imply that our solutions converge to a unique radial 



solution centered at a given point, since the drops may slowly move around a range of round profiles 
and may not converge to a single one. The physical uniqueness of the limiting profile and the character- 



ization of its center remain open at the moment (Though see Proposition 5.2 for the discussion on the 
uniqueness) . 
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2. The Equilibrium Problem 
Now consider the volume constrained minimization problem 

u = V}. 



inf{J(u) : M e iji(M^) and / u = V}. (2.1) 

Jm 

One can show immediately using symmetric decreasing rearrangements that any infimizcr of 

M{J{u) : u e H\R^) and / m = F and u radial } 

is also an infimizer of (2.1). That the only minimizers are radial is more delicate. One can show this 
using a theorem from [8]. We will not discuss the specifics here since this fact will not be needed for our 
arguments. 



To find the unique (up to translation) minimizers of (2.1 ) among radial functions first fix a particular 
radius r for the radially symmetric support set of the droplet 

{u > 0} = Br{0). 

Now the minimization becomes 

inf{ J(m) : u G HQ{Br{0)) and / u — V and u radial } 

the unique infimizer is just the solution of the Dirichlet problem for the Laplace operator with lagrange 
multiplier A = A(r) G (0, +oo) chosen so that the height profile 



Urix) = max{^(r2 



|xp),0} 



has the correct volume. Then explicitly calculating A(r) and minimizing J over r > one can easily 
check that for 

r ^ ( ^ )l/(N + l)yl/(N+l) 

* + 

Ur, is the strict minimizer of J among radial functions and therefore - by the rearrangement argument 
that was mentioned above - is also a minimizer among all functions. 

Alternatively, one can consider the Euler-Lagrange equation for (2.1) which is given by, 

{u > 0}, 



-Au = A 
A > 



s.t. 



(EQ) 



\Du\ = 1 on d{u > 0}. 
Then there is a classical theorem of Serrin [25^ regarding the uniqueness of solutions of (EQ). 



[0, +oo) be compactly supported such that d{u > 0} a 
hypersurface in M.^ and u is a classical solution of (EQ). Then {m > 0} = Br'{xo) for some Xq E M.^ 



Theorem 2.1. (Serrin) Let u 



and 



N + 2 



)l/(JV+l)^l/(W+l)^ 



(2.2) 



Serrin's proof of this theorem used the method of moving planes and a variant of the Hopf Lemma for 
domains with corners. The best possible spatial regularity we are able to show for the evolving contact 
line in (P-A| is Lipschitz. In order to apply Serrin's result to show the convergence to equilibrium, we 
need the same result to hold for viscosity solutions of {EQ) with positivity sets which are Lipschitz 
domains. It is not clear whether it is possible to use a variant of the moving planes method to show this. 



Instead we use some regularity results for free boundary problems. First a theorem of De Silva from 13 
shows that viscosity solutions of {EQ) with Lipschitz free boundaries are classical solutions. 



DYNAMIC STABILITY OF EQUILIBRIUM CAPILLARY DROPS 7 

Theorem 2.2. (De Silva) (Lipschitz implies C^'") Let C a domain and [0, +00) be a 

viscosity solution of the free boundary problem in f2, 

^ o-ijdijU — f in {x £ V, : u{x) > 0} 

(2.3) 

\Du\^g on r{u) := d{u > 0} n D, 

where Oij, g € C*''^(r2) for some < (3 < 1 and f G C{Q) H L°°(r2) and g > 0. Moreover, we have the 
ellipticity condition: there exists < A < A such that for all ^ € 9i?i(0) C and all x £ Vl, 

^<^a,j{x%ij < A. 

Then if Xq £ r(u), g(xo) > and T{u) is a Lipschitz graph in a neighborhood of Xq then Tiu) is C^'" in 
a smaller neighborhood of Xq for a > depending only on N , /3, and the ellipticity constants A and A. 



Using above result, higher regularity of u can be derived from the Hodograph method 22 when the 
coefficients are smooth and T{u) is locally Lipschitz. See Appendix C for more details . 



Corollary 2.3. Let u solve (2.3) with {uij) as identity matrix and f = g = 1. In addition suppose that 
T{u) is locally Lipschitz and is bounded. Then u and T{u) is C°° . 

3. Viscosity Solutions 



3.f. Basic Definitions and Assumptions. We will recall the basic viscosity solution theory for (P-A). 
A more detailed exposition for free boundary problems of a similar form can be found in 20 . First we 
restate the problem under consideration, 

-Au{x,t) ^ X{t) in fit, 

ut ^ Fi\Du\)\Du\ on Tf 

A : [0, +00) — > [0, +00) will be bounded and continuous. We will work in a space-time parabolic domain 
Q — U X {a,b] where < a < 5 < +00 and U is a domain (possibly unbounded) in M^, with parabolic 
boundary, 

dpQ ■.= U X {t^O}UdU X [a, b]. 
The positive phase of a height profile u : Q ^ [0, +00) is defined as 

n{u) {u > 0} and r(u) := dn{u) (3.f) 

and at particular times, 

fttiu) ^n{u)n{{y,s) £ Q : s ^t} and rt(u) = 9f2f (u). 

The dependence on u will be omitted when it is unambiguous which droplet profile we are referring to. 
The free boundary velocity F will satisfy the boundedness and ellipticity assumption, 

Assumption A. F : [0, +00) — > M is strictly monotone increasing and continuous, and F{1) = 0. 

Assumption B. There exists c > such that for e sufficiently small, 

{l + e)F{{l + ey^s) + ce>F{s). (3.2) 

Notice if Assumption 2 holds then necessarily F has sublinear growth at 00. For any s > using 
(3.2) [s/ej times with e — >■ yields: 

F{s) <cs + max{(l + s)F{0),0}. 



The monotonicity assumption on F implies, at least formally, that the problem (P-A I has a comparison 
principle. This is the underlying reason why viscosity solutions are the natural definition of weak solution 
for this PDE. The second assumption is important for proving the strong comparison results and thereby 
the regularity of the viscosity solutions of (P-A). ft is not clear to the authors whether this assumption 
is only technical. Some examples which satisfy the Assumptions \K\ and [B] are 
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Example 3.1. If F{s) := sP 
with c = 1. 



1 for p < 1 then by a simple calculation (3.2) is satisfied for all e > 



The following example in the case p = 2 will be important to us later, 



Example 3.2. Let M > and p > 1, suppose F{s) := maxja^* — 1,M} then (3.2 1 is satisfied for 
e < 1 + i(p — I)~^ and c > 2(M + 2). The calculation is given below, 

max{(l + e)[(l + e)-PsP - 1)], M} + ce > Taax{sP - I - (p - l)esP - e, M} + ce 

> max {sP -1- 2{M + 2)e, M} + ce 

> max{sP - I + (c - 2(Af + 2)) e, M} 

> max{sP - 1,M}. 



Above in the second line we have used that if 



I 



then 



- (p - l)esP -€<M 
sP < 2{M + 2). 



Now we turn to defining a notion of solution for (P-A). For E C x [0, +oo) we use the notation 
C'^'^{E) for functions / with two continuous derivatives in the spatial variables and one continuous 
derivative in time. First we define a classical solution of the free boundary problem. 

Definition 3.1. (Classical Solutions) A profile m : Q — ?> [0, +oo) is a classical subsolution (supersolution) 
for the free boundary problem (P-A) on Q if m G C^'^(il(M) n Q) and (P-A) is satisfied in the pointwise 
classical sense with equality signs replaced by < (> in the supersolution case). Classical solutions are 
both super and subsolutions. 

For general initial data the contact line motion problem will not have a classical solution. To get 
existence of a weak solution we define the notion of viscosity solution of ( P-A ) . First though we define 
the following standard notion. 

Definition 3.2. (Strictly separated) Let u, w be defined on a set _D C then we say v and w are 
strictly separated on D and write t; -< w if D,{v) H 13 is compact and v < w on Q.{v) H D. 

Next we define subsolutions and supersolutions, then a viscosity solution is defined to satisfy both 
the sub and supersolution properties. Informally, m is a subsolution of (P-A ) if u cannot be crossed from 
above by any strict classical supersolution. More precisely: 

Definition 3.3. (Subsolution) A non-negative upper semi-continuous function u : Q ^ M+ is a subso- 
lution of (P-A) if, for any parabolic neighborhood Q' C Q, and any strict classical supersolution </) with 
u < 4> on dpQ', then u < (0)+ in Q' . 

Definition 3.4. (Supersolution) A non-negative lower semi-continuous function u : Q ^ M+ is a 
supersolution of ( P-A[ ) if, for any parabolic neighborhood Q' C Q, and any strict classical supersolution 
(j) with u ^ (j) on dpQ' , then (j) < u in Q' . 



Definition 3.5. (Solution) A viscosity solution of (P-A) is a non-negative continuous function u on Q 
which is both a supersolution and a subsolution. 

Naturally, one can assign boundary data on the parabolic boundary. We will usually have Q = 



pN 



(0, T] and in that case this will reduce to assigning initial data. 



Definition 3.6. (Subsolution with boundary data g) Let g : dpQ — ?► [0, -l-oo) bounded then u is a 
subsolution of ( P-A| ) on Q with boundary data 5 if u is a subsolution of (P-A) on Q as defined above 
and 

lim sup u{y, s) < g{x, t) 
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Supersolutions and solutions are then defined analogously. 

For the rest of the paper, in the case Q = x (0, T], when we assign initial data we will specify the 
positivity set as some J7o C an open domain and then ug the initial data will be taken to solve 

— Auo(a;) = A(0) and uq{x) = for a; G Fq = Silp- 

As mentioned in the introduction, it is often more natural to think of the problem as an evolution on 
the positivity sets. 



Remark 3.3. In general, the solution of (P-A I is not expected to be continuous. First of all, the solutions 
can vanish or blow up in finite time. The following example in the case = 1 demonstrates these 
behaviors. 



Example 3.4. Consider in iV = 1 the problem (P-A) with X{t) = A for some A > and initial data 

wo(a;) = 2^(l-x2) + , r!o = [-l,l]. 



The solution of ( P-A ) with the above initial data and inhomogeneity will take the form 

u(x,t) = ^[r{t)~x^)+ 
where r[t) solves the ordinary differential equation 

f{t)^F{^r{t)) 

with initial data r(0) = 1. From basic ODE theory there are choices of F for which some solutions of 
the above equation blow up in finite time. For example if 

F{s) = s2 _ 2 and A > 2N, 

then first of all r{t) > as long as the solution exists and moreover 

f(t) > Ar{tf - 1 > 3r(i)2 

an equation for which finite time blow up is well known. (This behavior is ruled out by Assumption 2 on 
F which implies sublinear growth at oo.) Alternatively, even when F does not have super linear growth 
solutions can disappear in finite time. Consider now any F such that 

F(0) < 

then let —F{{)) > (5 > and choose A sufficiently small that 

F{^)<-5<Q. 
Then we have the following differential inequality for r(i), 

f{t) < ~S 

so the solution u{x,t) must become extinct before time 1/5. 

Even without blow up the solution may be discontinuous for certain initial data: one can construct 
an example in one dimension where two adjacent drops merge after a short time, causing a discontinuity. 
To address the long time behavior of our solutions without the complexity of multiple components, we 
will make restrictions on the initial data such that our solutions actually are continuous. 

The following lemma clarifying the connection between classical and viscosity solutions is standard 
in the viscosity solutions theory. For example see [6 20 . 



Lemma 3.5. Suppose that u : Q ^ [0, -l-oo) is a classical suhsolution (supersolution) of (P-A) in Q 
then it is also a viscosity suhsolution (supersolution) in Q. Conversely suppose that u is a viscosity 
suhsolution (supersolution) of (P-A) in Q and moreover u is sufficiently regular, in particular 

then u is a classical suhsolution (supersolution) in Q. 
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An important basic property of viscosity solutions theory is the stabihty of viscosity solutions under 
uniform convergence. We state this in the following lemma. For example see [?]. 



Lemma 3.6. Let u„ 



^ X [0, T] be a sequence of viscosity solutions of the problems 
Au„{x,t) ^ Xn{t) in ^t{Un), 

FJlDUr^ 



dtUn ^ Fn{\DUn\)\DUn\ On Tt{u„). 

where Fn are all monotone increasing and continuous. Suppose that u- 
uniformly on compact sets, then u is a viscosity solution of 

—Au{x,t) — X{t) in flt{u), 
dtu = F{\Du\)\Du\ on Tt{u). 



u, Fn 



F and X„ 



As mentioned in the introduction, it is useful sometimes to formulate the problem (P-A ) as an evolution 
. We quickly define what it means for an evolution 

Qt ■■ [0, +oo) { bounded subsets of M^} 



on domains in 



to solve (P-A I in the viscosity sense. 



Definition 3.7. The evolution fit is a supersolution (subsolution) of (P-A) if and only if, 

-Au{x,t) = X{t) in fit 



(3.3) 







is a supersolution (subsolution) in the sense of Definition |3.4| (Definition 3.3 1 



3.2. Sup and Inf Convolutions. An important property of subsolutions (supersolutions) is the closure 
under sup (inf) convolutions. These convolutions will be used to perturb the originial solutions and to 
overcome the lack of scaling invariance in the problem ( P-A I . 



Lemma 3.7. (Convolutions in space) 

(a) Let u be a subsolution of (P-A I on . 
ofu 



dN 



X [Oj-l-oo), and r > 0, c > and define the sup-convolution 



u{x,t) := sup u{y,t). 



(3.4) 



Then u is a subsolution with free boundary speed F(\Du\) — c as long as r ^ ct > 0. 
(b ) Let u be a supersolution of ( P-A I on X [0, -|-c»), and r, c > and define the inf- convolution of u 



u{x,t) := inf u{y,t). 



(3.5) 



Then u is a supersolution with free boundary speed F{\Du\) -\- c as long as r — ct > 0. 

Before we prove this let us note the continuity properties of the sup and inf convolutions. 

Lemma 3.8. (Continuity of sup and inf convolutions) Let f : ^ R and let 5 : ^ B{R^) such 
that S{x) is compact for all X eR^ and S is continuous under the Hausdorff distance dn ■ Define the 
sup and inf convolutions 

f^{x)^ sup f{y + h) and fs{x)= ini f{y + h). 
hes{x) hes{x} 

Then we have the following 

(i) If f is upper semi-continuous then f^ is as well, 
(ii) If f is lower semi- continuous the fs is as well. 
(Hi) If f is continuous then both f^ and fs are as well. 
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Proof. We will prove (i), (ii) is similar and (iii) follows from (i) and (ii) by noting that suprema of 
lower semi-continuous functions are lower semicontinuous and infima of upper semicontinuous functions 
are upper semi-continuous. Let a € M we will show that the sub level sets {f^ < a} are open. Let 
xq £ {f^ < a} then xq + S{xo) C {/ < a}. Since xq + S{xo) is compact, {/ > a} is closed, and they 
are disjoint, they must be a positive distance apart. That is, there exists e > such that 

U B,ixo + h)C{f <a}. 

heS(xo) 

By continuity of S there exists S < e/2 such that \x — xo\ < S means that dH{S{x), S{xo)) < e/2. Let 
X e i?£/2(a^o) a-nd since upper semi-continuous functions achieve their maximum on compact sets there 
exists h S S{x) such that {x) — f{x + h). Then there exists ft-o in *S'(xo) such that \h — /iqI < e/2 so 
that 

|(a; + h) — (xo + /io)| < e and hence x + h £ B^{xo + ho) C {/ < a}. 
So f^{x) < a and therefore B^f2{xo) ^ {f^ < ct}. 

□ 



The proof of Lemma 3.7' We will only prove the inf convolution case (b), since the proof of (a) is 
essentially the same. As in (3.5) for r, c > define the inf convolution of a supersolution w, 

u(x, t) := inf u{y, t). 

Suppose that for some parabolic domain Q C x [0, r/c) there exists a strict classical subsolution 
(p € C^'^{Q n with the free boundary speed F{\D(I)\) + c, which satisfies 

(j) ^ u on dpQ and (f){xo,to) > u{xo,to) for some {xo,to) e Q. 

Let xi be the point where u(xi,to) = u{xo,tQ). Now we define the translated (and convoluted) test 
function 

■0(x,t):= sup (f){y + {xq — xi),t) tor t < to, 

\x-y\<c{to-t),(y,t)eQ 

and define 'ip2 by the solution of the Dirichlet problem 

-A^2i-,t) = \{t) in > 0} n Q, ■02 = -0 on the lateral boundary of {■(/;> 0} n Q. 

Note that the free boundary speed of ip2 has been decreased by c over the free boundary speed of (j) in 
Q. In other words, ^p2 is a classical subsolution with free boundary speed F{\Dijj\). 

Due to the definition of u, we have 02 ^ on the parabolic neighborhood of Q -I- (xi — xq) and "02 > " 
at {xi,to), which yields a contradiction. 

□ 

3.3. Comparison. We state the general strictly separated comparison principle for viscosity solutions 
of (??), 

Theorem 3.9. (Comparison for strictly separated data) Suppose u is a supersolution and v a subsolution 
of (P-A) on Q. Suppose u and v are strictly separated (defined below) on the parabolic boundary of Q, 

V ^ u on dpQ 



then V < u in Q{v) D Q, in particular T{v) cannot touch T{u) from the interior on any compact subset 

ofQ. 

Proof. The proof can be found in various places in the literature [6j[l0j|20|[2l] . We will go into more 
detail later in Section l3?3l □ 



As mentioned in Remark |3.3[ we will need to make a geometric restriction on our initial data to expect 
the existence of a viscosity solution which is stable under a family of perturbations. First, we recall the 
definition of a set star-shaped with respect to a point: 
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Definition 3.8. A domain fl C M. is called star-shaped with respect to a point x if for every y g dil 
the line segment between x and y is contained in 

Then the assumption on our initial data fig will be called strong star-shapedness: 

Definition 3.9. We call a domain 57 strongly star-shaped if there exists r > such that ft is star-shaped 
with respect to every point of Br{0). For each i? > we define the class of uniformly bounded strongly 
star-shaped sets: 

Sr,B. ■= C Bfi{0) : V, star-shaped with respect to Br{0)}, S,. := Ui?,>o '^'-.fl' (3-6) 
and the class of strongly star-shaped sets 

So := U Sr. (3.7) 

That the ball -6^(0) is centered at the origin is only for convenience, the problem is translation 
invariant. We will often refer to the strongly star-shaped property by saying that a set fig € Sr,R in 
order to clarify the role of r. We note a basic property of sets in Sr,R ~ they are Lipschitz domains with 
the Lipschitz constant depending only on r, R. First, we define a notation for cones with apex at the 
origin. For X e and 61 e (0, vr) the cone in the direction x with opening angle 9 is called, 

C{x,0) ■.= {y:{x,y) > (cos 9)\x\\y\} . (3.8) 

We show that Cz Sr have interior and exterior cones at every boundary point: 

Lemma 3.10. The following are equivalent for a bounded domain D, D Br{0): 
(i) The domain J7 € 5^. 

(ii) There is an Cq G (0,oo] such that for all x G ft'~' there is an exterior cone to il, 

x + C{x,0x)r\B^g{x) cfl^ where sin6'^ = (3.9) 
(Hi) For all x e dfl there is an interior 'cone' to 

{x + C{-x,9J)nC{x,^-0^)UBr{O)cn where sm9^ ^ (3.10) 
(iv) There exists Eq > so that 

f7cc Pi [{l + e)n + z] for all eo > e > 0, (3.11) 

\z\<ae 

for every < a < r. 

Proof. The proof for parts (i)-(iii) is essentially in the picture, see Figure [2j For part (iv) let us first 
note that (see figure [2|, 

X + C{x,0^)nB,„r{x) C y U Baeiil + e)x). 

£o>e>0 0<a<r 

Therefore by (ii) 57 e 5^ if and only if for all x £ ft'-'' all eq > e > and all < a < r, 

Baei{l + e)x) CC n^, 
or equivalently, for every cq > e > and all < a < r, 

{x : d{x, (1 + 6)57*^) < ae} CC . 
Again, equivalently, for every eq > e > and all < a < r, 

ncc{x:d{x,{l + e)n^)<ae}^=R^\ [j [{1 + e)n'^ + z] ^ f] [{I + e)n + z] . 

\z\<.ae |z|<ae 

□ 
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Figure 2. Interior and exterior cones for a strongly star-shaped domain. 

We prove a comparison principle for reflections through hyperplanes. This will be very useful to us 
later because the reflection ordering is a property which does not depend on the radius from the strong 
star-shapedness property. Let H he a hyperplane in with y £ H define the reflection through H by 

(I)h{x) =x-2(x-y, vh)vh- 

The symmetry of the problem ( P-A[ ) with respect to reflections at least formally implies that if a solution 



u{-,t) and its reflection u{(j)H{'), t) through H are initially ordered in the half spaces of H then they will 
remain so. 



Proposition 3.11. (Reflection Comparison) Suppose [0, T] — > [0, +oo) is a solution of ( |P-A[ ) 

such that flt{u) € Sr for some r > and for all t G [0,T]. Let H he a hyperplane in with 
H n Bj,{0) = and define and H_ the half spaces of H such that i?r(0) C Suppose that 

(t>H{^o) n i/+ C fio n i/+ i.e. 

u{4>h{x), 0) < u{x, 0) for X e H+. 

Then we have the ordering: 

u{(f>H{x),t) < u{x,t) for {x,t) e H+x [0,T]. 
A similar argument will give the more standard strong comparison result, 



Lemma 3.12. (Strong comparison) Suppose u and v are respectively a sub- and supersolution of (P-A| 
on X [0, +00) which are initially ordered, 

no{u) c no{v). 

Suppose that that there exists T > such that for each t € [0,T] either flt{u) or Qt{v) is in Sr^R for 
some r,R> independent oft. Then let c from Assumptions!^ and define to := r'^/{Rc) then 

nt{u) C nt{v) fort<T + to. 



Proof. The proof is a slightly easier version of the proof of Proposition 3.11 so we omit it. □ 



In order to get a strong comparison type result we want to slightly perturb the supersolution (or the 



subsolution) so that strict separation holds initially and we can use Theorem 3.9 To achieve this we 
use the inf and sup convolutions introduced above. This will be where the 'technical' assumption on F 
iBl comes in. 
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Proof of Proposition \ 3. 1 i} ' Let eq ^^nd c from Assumption [Bj < a < r, and define for < e < eo, 

We(x, t) = (1 + e)^u((l + e)"^a;, i); and 
Ue(x,t) = inf v^{y,t) for < i < -. 



Then, from Lemma 3.10 part (iv) and the strong maximum principle for harmonic functions, we have 

u(x,t) < uJx,t) in X [0. -) for e > 0. 

c 

We check the super-sohition property of u^. Note that Vc{x,t) is a supersolution of, 

-Av^{x,t) > X(t) in nt{v,) 

dtV,>il + e)F{{l + e)-'\Dv,\)\Dv,\ on Tt{v,). 
Then, from Lemma 3.7 Ue(-,t) is a supersolution of 

-Au^{x,t) > X{t) in fJt(uc) 

dtu,>[{l + e)F{il + e)-^\Du,\)+ce]\Du,\ on rt(v,). 



Using Assumption [Bl we conclude that is a super-solution of ( P-A I . 



Since is a super-solution of (P-A), by the strictly separated comparison principle Theorem 3.9 



u{(j)H{x),t) < u^{x,t) in i/+x[0, -). 

c 



Since 



u as e — >■ 0, it follows that 



u{(l)Hix),t) < u{x,t) in _ff+x[0,-). 



We now iterate [Toc/a\ + 1 times to conclude. 



□ 



3.4. Short Time Existence. As usual in the viscosity solutions theory, the comparison theorem is the 
key ingredient needed to use Perron's method to show existence. In the existence proof we will use the 
following elementary facts about solutions of Poisson's equation at the boundary. 

Lemma 3.13. (Boundary Holder Estimates in Lipschitz Domains) Let ft C M.^ be a domain and 
xq £ dil and suppose that there exists 9 e (0,7r/2) and v G S^^^ such that, 

xo + {x: {x, z/) > cos 6l|a;|} C \ n. 

In other words, Q has an exterior cone of opening angle 6 at xq. Suppose A > A > and u : H — )• [0, +oo) 
satisfies 

—Au{x) < A for X Cz Q 

u(x) — on dfl 

Then there exists a — a{6, N) e (0, 1) and C ~ C(A, N) > such that for all h > 0, 

sup u{xo + z) < eh"'. 

\z\<h 

Lemma 3.14. (Boundary Gradient Estimate in C^'^ Domains) Let C be a domain xq € d^l and 
suppose that there exists r > 0, y € M.\n such that 

Br{y)nU = {x}. 

In other words, fl has an exterior ball at Xg of radius r. Suppose A > A > and u : f2 — > [0, +oo) 
satisfies u{xo) = and 

—Au < A for X Cz il. 
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Then 



lim 



sup 



We move on to the 'short time' existence theorem. Actually this shows the existence of a global in 
time discontinuous viscosity solution without any need for Assumption[B] The uniqueness and continuity 
however rely on the strong comparison result which requires Assumption [B] and only holds for a short 
time. This is to be expected as noted in Remark 1 3. 3 1 In particular the limiting factor will be the strong 
comparison principle Lemma [3 . 1 2| which only holds for a short time barring any a priori knowledge about 
strong star-shapedness. 

Theorem 3.15. (Short time existence and uniqueness) Let r, R > and CIq G Sr,R and let uq solve 

—Auq{x) = A(0) and Uo{x) — for x G Tq. 

Then there is a to > depending only on r such that there exists a (unique) continuous viscosity solution 
u of (P-Al on Q — x (0,io] with initial data uq. 

Proof. This proof can be found in 14' . We construct a sub and supersolution which take the initial data 
continuously then apply Perron's method. 

1. First we construct a subsolution. Let 7 > to be chosen later and define, 

(1-70'«o(t^) for i<7"^ 



U{x,t) 

for t > 7~^ 

Then, calculating formally, if xq G Fq then (1 — jt)x() is in Tt and 

Ut ' 



(3.12) 



r((l - -ft)xo) = 7X0 



if 7 = "f{r,R) is chosen small. In particular, 



<-7|a;oK 1- 



r 
\xo\ 



< mini^ 



7 < mini^ (R^ - r^) 



-1/2 



will work. This calculation shows that J7 is a subsolution of (P-A) when dflo is smooth enough that 
Duq is defined on dflQ. The calculation can be transferred to the test functions to show that J7 is a 
subsolution in the general case. 

2. Next we construct a supersolution which takes the initial data continuously. Let a = a{r/2,R) G 
(0, 1) from Lemma 3.13 combined with Lemma 3.10 Define, for any h > 0, 



n ii + h)no 

\z\<p{t) 



where 



Vhix,t) = u 



p{t) =a[h~ - max{F (C/i""^) , l}n for < t < max{F (C/i"^^) , 1}-^ 



(3.13) 



Here < a < r so that mo -< Vh as in Lemma [3.10| It follows from Lemma [XT] that Vh satisfies 

dtVh 



\DVj, 



-{x,t) > max{F {Ch°'-^) , 1} on Tt{Vh) 



in the viscosity sense. Let x G Ff(V;i), then there exists y G dBp(t){x) H (1 + h)To and ilt{Vh) has the 
exterior ball at x 



Then, for h sufficiently small, ^lt{Vh) is an intersection of domains star-shaped with respect to 5^/2(0) 
and we apply Lemma [3. 14| and then Lemma [3.13| to get, 

\DVh\ix,t) <Ch-^ sup Vh{x + z,t) < Ch'^-^ 

\z\<2h 
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and it follows from the monotonicity of F that Vh is a supersolution. 
3. Now we apply Perron's method. Let us define 

u{x,t) := sup{w : x [0, +00) — > [0, +00) : u is a subsolution of (P-A) with initial data uo}- 

From its definition and Theorem |3.9| we get the ordering for all /i > 

U{x,t) <u{x,t) <Vh{x,t) for < i < max{F(C/i"-i) (3.14) 

Let us now define the upper and lower semi-continuous envelopes of u, respectively 

u*{x,t) = lim sup u{y,s) and u■^,{x,t) — lim inf u(y,s). 

It is standard in viscosity solution theory (see [12| ) to show that u* is a subsolution and is a super- 
solution. Due to (3.14) one can check that u{-,t) — >■ uq{-) uniformly as t — > 0, i.e. u* = at t = 0. 
Thus, letting to — to^r, R) > as in Lemma 3.12 we have u* < m* on [0,to]- Combining this with the 
fact < w*, we obtain the existence of a unique continuous viscosity solution of (P-A I on [0,fo]- 

□ 

We point out that the above proof yields a modulus of continuity in time of the contact line at t = 0. 
For solutions which are uniformly strongly star-shaped on a time interval [0, T] this extends to give a 
uniform modulus of continuity in time for the contact line on [0,T]. We make this explicit in the case 
when F has polynomial growth below. For the rest of this section we will make the following assumption: 



Assumption C. There exists p > such that: 



lim sup < +00 



(3.15) 



We prove a modulus continuity for strongly star-shaped solutions of (P-A) depending only on r, R, 
p, sup A and s\vpF{s)/s^: 

Corollary 3.16. Suppose [0,r] — > [0, +00) is a viscosity solution of ( |P-A[ ) with < \(t) < A 

bounded and F satisfying Assumption s\A\ \b\ and^^ Moreover suppose that f2t(w) G Sr.R for some r, R > 
and all t G [0,T]. Let a from Lemma 3.14 and 

^ ■ 1+J9(1-q) ' 

then Tt{u) is in time in the H aus dor ff metric, i.e., 

dH{TtAu).^t,{u)) < C2\ti-t2f for alltiM G [0,r]. 
Here C2 depends only on svLpF{s)/.sP, p, A and N. 

Remark 3.17. An analogous equicontinuity result is true for general F satisfying Assumptions \K\ and 
[B) In this case the Holder modulus of continuity will have to be replaced by a more general modulus 
of continuity depending on the growth of F at infinity. We restrict to the polynomial growth case for 
simplicity of presentation. 



Proof. Let T > i2 > > and 7 from the proof of Theorem 3.15 Using the barriers from the previous 
theorem, 

(l-70VT^,ii) 



for i — <i < 7 
for t — ti > 7" 



-7* ' 

u{x,t2) > U{x,t2) := 

Also for h = C(<2 — ti)^ we have that 

uix,t2) < Vh{x,t2) < il + h)^u 

Then the following containments hold: 

[(1 - 7(i2 - <i)) V 0]^! (u) C a. (u) C (1 + C{t2 - tiY)nt, {u). 



(3.16) 



(3.17) 
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So if a; e Ff ^ (u) and (^2 — ^i) small, 

{y-\y-x\<CR{t2~hf}nTtJu)^' 



and 



dH{rtAu),rt^{u))<c\h~hf. 



(3.18) 

□ 



Remark 3.18. See Lemma A. 3 for some implications of (3.18). 



We will be able to use this Corollary to show existence of a continuous viscosity solution for (P-Al 
when F has polynomial growth. The idea is to take the unique solution um for the problem with the 
free boundary velocity max{i^, M} and let M — >■ oo. Corollary 3.16 will give equicontinuity to take a 
convergent subsequence of the um as long as the fltiuAi) € Sr,B. with r, R uniform in M and t. In order 
to show this we will need some kind of preservation of the strongly star-shaped property which does not 
depend on M. 

Lemma 3.19. Let T > and u„ : M.^ x [0,T] — > [0, +oo) be viscosity solutions of: 

~AUn{x,t) ^ Xnit) in fit(Un), 

dtUn = Fr,{\Dun\)\Dun\ on TtK). (3.19) 

Here A„ — > A > uniformly so that in particular there exists A > so that: 

supsupA„(t) < A < +0O. 

n t 

We suppose that < F for all n and F^ ~^ F uniformly on compact sets. The limiting free boundary 
speed F has polynomial growth of order p > 1, 

F{s) 

< lim sup < +00. 

s-^oo S^ 

The Fn are assumed to all satisfy Assumptions^^^^and^^ with Assumption^^ satisfied uniformly : 

sup„ sup^ < +00. (3.20) 

Moreover, suppose that there exist r,R> such that for all n and all t e [0,T] we have Qtiun) G Sr.R- 
However note that we do NOT assume that the constants c„ from Assumption\B\ are uniformly bounded. 
Then the w„ converge uniformly on X [0, T] to a viscosity solution of: 



-Au{x,t) ^ \{t) 
dtu = F{\Du\)\Du\ 



in 51* (u), 
on Tt{u). 



(3.21) 



// additionally A„(<) < X(t) for all n then u is the minimal viscosity solution (3.21). In particular, in 
this case, the limit does not depend on the approximating sequence Fn. 



Proof. 1. First we show that f2t(u„) is an equicontinuous sequence of paths in Sr^R- From Lemma A. 3 
and Corollary |3.16| we derive for i, s > 0: 

S\XpdH{^t{Un),^s{Un)) < sup d// (Fj (u„) , F^ (u„) ) <C\t-sf. 



Then from the compactness lemma A. 4 for paths in Sr.R up to a subsequence VLt{un) converge uniformly 
to a continuous path VLt in Sr,R, with free boundary F^ :— dflf. 
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Let u{x,t) be the solution of 



-Au{x,t) = A(t) for a; e fit 
u{x,t) = for x^Tt- 



(3.22) 



We claim that u„ converges uniformly to u in x [0,r], we will demonstrate this by showing the 
following estimate: 



u\{x,t) 



< 



r,R 



dH{ntM,nt{u)r + \Xr.{t)-x{t)\. 



Let X e r2t(u)Af2f (u„), then: 

d{x, dintiu) U 17t(u„))) < dniTtiun), Tt{u)) S,fl dif(r!t(M„), r2t(u)). 
Since m„ — u = on d{il,t{u) U ilf(u„)) we get from Lemma [3.13[ 

\Un - U\{x,t) < Cd{x,d{^lt{u) U ^t{UnW <r..R dH{^t{Ur,) Mt{u)T ■ 

Meanwhile for x G ^lt{u) H rit(it„) we combine the above inequality which holds on the boundary of 
n Vlt{un) with the fact that diam(17t(u) n f2i(u„)) 1 to get, 

\Un - U\{X, t) < dH{nt{u„),ntiu)r + \Xnit) - X{t)\. 



We apply the stability of viscosity solutions under uniform convergence, Lemma 3.6 to see that u is a 



viscosity solution of the PDE (3.21 1 as claimed. 



2. Now we show that if A„ A then u must be the smallest viscosity solution of (3.21 1. Let v be 



another viscosity solution of (3.21 1. Note that due to the orderings F„ < F and A„ < A we have that v 



is a supersolution of each of the approximating problems (3.191. Using the strong comparison principle 



Lemma 3.12 which holds for each problem (3.191 we get 



Un < V on 



[o,r]. 



□ 



3.5. Preservation of the Strongly Star-shaped Property. Here we describe some of the properties 



of the viscosity solutions of (P-A) and (P-V). First, we will show that, if the droplet initially is in Sq, 



then this property persists for a short time with the radius of the strong star-shaped property going to 
zero in some finite amount of time. This short-time regularity will not be sufficient to prove any kind of 
long-time behavior. However, in Proposition |3. 1 1] we have shown that a reflection comparison principle 
holds as long as the positivity set is in Sr^R for any r, R > 0. The reflection comparison principle in 
turn will allow us to show, in some situations, that actually there was no loss of star-shapedness. In 



particular for ( P-V I we will show that initial data which is in Sr for a sufficiently large r along with a 



condition on the smallness of 



sup \x\ 



inf \x\ 



will maintain some regularity globally in time. In fact, it will be strongly star-shaped for all time with 
a possibly smaller radius. 

We show that strong star-shapedness cannot disappear immediately. The below Lemma is essentially 
contained in jl4j, we use a different proof. 

Lemma 3.20. (Short time strong star-shapedness) Let [0, -t-cxi) — > [0, +(X)) be the solution of 

(P-A) with initial positive phase Oq G Sr for some r > 0. Then 



r-ct for 0<t<- 
c 



where c is from Assumption\E[ 
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Proof. Let eg > e > and Q < a < r, then from Lemma |3. 10| and Assumption [B] we know that 



u,{x,t) ■.^Mi,i<ae-cet{l + efu{{l + e)-\x + z),t) for i<f 

is a supersolution of ( |P-A ) which has u(-,0) -< u^{-,0). Therefore from the strict comparison result 
Theorem 3.9 u < for t < a/c and in particular, 

Ot(u) C Ot(ue) = Pi ((1 + e)rit(M) + 2:) for every eo > e > 0. 

|^|<(a — ci)€ 



Since e and < a < r were arbitrary the converse direction of Lemma 3.10 part (iv) implies that rit(it) 
has the claimed star-shapedness. 

□ 

We would like to show ilt {u) G Sr as long as Br (0) is contained in the positive phase Vtt (u) . This kind 
of preservation of the star-shapedness property would be very useful because we could get regularity 
results by simply showing that the positivity set must always contain some small ball around the origin. 
However, the arguments of Lemma 3.20| are insufficient to prove such a result. This is where the reflection 
comparison principle Proposition |3.11 is useful. The key fact is that reflection comparison holds as long 
as r2t(u) G Sr for any r > 0. 

We begin by defining a slightly stronger property than strong star-shapedness which is defined in 
terms of refiections. In order to do this we will need some notations. Let i7 be an open bounded domain 
in M^. For 1/ G S^^^, let H — H{v) be the hyperplane through the origin orthogonal to v with the half 
spaces it defines: 

H+ = H+{v) =^{x:x-v>0}., iJ_ = H^{v) = {x -.x-v <0}. 
Now for s > define the translates, 

H{s) + siy, H+{s) = H+ + sv, (s) = + sv. 

For s sufficiently large 

VL C i?_(s) 

and so trivially 

(l)H{s){^r\H+{s)) (znr\H^{s). (3.23) 

See Figure [3] for a depiction of the situation. Now let us slide the hyperplane inwards towards the origin 
by decreasing s until (3.23) no longer holds. We will call s,„in(:/, Vl) the closest to the origin that we can 
move this plane with (3.23) always holding. More precisely, we define 

Smin{v, n) := inf{s > : (t)H(t){^ H H+{t)) C n n H^{t) for all t > s}. (3.24) 

In the following we omit the dependence of Smin on when it will not present any confusion. If 
Smin('^) = for cvcry direction v then must be a ball. The basic approach of Serrin in j25) to showing 
the symmetry of solutions of (EQ) is to show that Synin{i^) — for all ly. In our case it is useful to use 
this same idea of symmetry but in a weaker form. 

Definition 3.10. We say that a bounded, open set ^ has p-reflection if -Bp(O) C and 

sup Smin(t^) < P- (3.25) 

Remark 3.21. As with strong star-shapedness, given some initial data flo, we will in general assume that 

sup Smin(i^, ^^o) = mf SUp S^in (l', ^^0 + ^) • 

Of course if this is not the case it can be corrected by a spatial translation. 
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Figure 3. Reflection comparison will hold. 



Remark 3.22. One can show that (see Lemma A.ll, upon fixing a maximal diameter i?, a set has p- 
reflection if its boundary normals deviate from the radial direction by 0{p/ R) and additionally the 
oscillation of the boundary is 0{p^/R). Although we would like to know whether the flow preserves 
convexity, the p-reflection property does not seem to be very useful in this direction. In particular 
/9-reflection does not imply convexity, and neither does the converse hold. 



Note that the property (3.25) is preserved over time due to the reflection comparison, as long as the 
positivity set ^lt{u) £ Sr for any r > 0. We list some basic facts about sets which have p-reflection. The 
proofs are postponed to the appendix. The first says that the condition p-reflection imposes a condition 
on the spatial location of the boundary of the set. 



Lemma 3.23. Suppose has p-reflection, then 



sup 

xedn 



inf \x\ < 4p. 

xedn 



(3.26) 



We mention that (3.25 1 imposes a bound on what kind of normal vectors dn can have. In fact D, will 



be strongly star-shaped with radius depending only on dist(9f2, -6^(0)). 

Lemma 3.24. Suppose f2 has p-reflection. Then fl satisfies the following: 
(a) for all x e 90 there is an exterior cone to f2 at x, 

x -t- C (x, 02,) C \ where cos(/)3. = A, e (0, 7r/2), 



(3.27) 



and C{x,(l)x) is the cone in direction x of opening angle (px as defined in (3.8); 
(b) U, ^ Sj. where 



r = r{p, mi \x\) = ( inf |a;|^ 



Now, combining the strong star-shapedness from Lemma [3.24| with the reflection comparison principle 
we obtain that p-reflection is preserved as long as the -Bp(O) is contained in evolving 



3.11 



Proposition 
positive phase. 

Lemma 3.25. Suppose u 
some p > 0. Let I = [0,r) he the maximal time interval containing on which Bp(0) C flt{u). Then 
r2f(u) has p-reflection and in particular flt{u) € Sr for some r > for every t £ I. 



oN 



X [0, +oo) — )• [0, +oo) is a solution of (P-Al and flo has p-reflection for 
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Proof. Note that by comparing with a radial subsolution placed below m(-,0) one can see that T > 0. 
Suppose towards contradiction that the lemma fails, that is 

t^, := inf{i G [0,T) : ilt{u) does not have p-reflection } < T. 



Since Bp{0) C ilt, (u) we have that 



h dist{nt,{uf,Bp{0)) > 



so that by Lemma 3.24 f2t^ (m) is star-shaped with respect to S(;j2^2p/i)i/2 (0). Then by Lemma 



3.20 



fit(u) e Sr for r = + 2pft.) 



1/2 



(3.28) 



on some slightly larger interval [t^,t^+r/c). Now let iJ be a hyperplane which does not intersect i?p(0). 
We can apply the reflection comparison Proposition 3.11 to u{x,t) on [i*,<* + r/c) to see that, 



u{(t)H{x),t) < u{x,t) for {x,t) e H+n [<*,i* +r/c) 

and therefore 

This holds for all admissible H, so r2t(u) has p-reflection on +r/c), contradicting the definition of 

U. □ 



Notice that an immediate consequence of Lemma |3.25| is that for any initial data which has 
p-reflection the evolution fit has p-reflection for at least as long as, 



t < 



min^gro F| - p 
mini^ 



The time on the right hand side above is independent of c from Assumption [b] and of X{t). 

Now we present an application of this idea to the volume preserving problem, showing that for certain 
initial data any solution of the volume preserving problem must have p-reflection for all time. One can 



think of this as an a priori estimate for solutions of (P-V). If the initial data flo has p-reflection then 



Lemma |3.25| says that the solution fit of the volume preserving flow also has p-reflection until such a 
time as touches Bp{0) from the exterior. For a domain il C let us define the associated Lagrange 
multiplier. 



X[n] ^ inf{ J \Dv\^ : v G H^{n) and J v = V}. 



Then from Lemma 3.23 at the touching time we know that fit, C B^p(0). This allows us to get a lower 
bound on the Lagrange multiplier A [fit J and show that Bp{0) is a strict subsolution near the touching 
time in the sense that 

M(P^-N^) 

2N ^ I I ; 

is a strict classical subsolution for \t—t.^, \ small. Of course, this will be a contradiction of the supersolution 
property of fit . We make this precise in the following Lemma. 



0,T] — > [0, -|-oo) is a solution of (P-V) with initial data that has 



Lemma 3.26. Suppose u 
p-reflection with Bp <Z<Z fl{). Then there is a dimension constant Cn such that if 

< p < CnVt^ (3.29) 
then there exists a > such that i?(i_|_a)p(0) C fltiu) for all t > 0. In particular, we establish that 

1 /N{N + 2y 



Cn > 
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Remark 3.27. The scaUng p ^ 1^"+^ is the natural one for such an inequality. In particular, if p is larger 
than the unique radial stationary solution, 



N 



N + 2 



V "+i 



then we can construct a counter-example using radially symmetric solutions of (P-Vl. The solution of 
(P-V) with initial data Br{0) for any r > p has p-reflection and converges as i — >■ oo to Br,{0) so, in 
particular, there is some time after which i?p(0) is no longer contained in the evolving positivity set. Of 
course in this case the initial data has p-reflection for every p > 0. 



Proof. From the compact containment Bp CC ilo along with the strict inequality (3.291 let a > small 
so that: 



S(i+a)p(0) CC ilo and p < 



N{N+2) 
(l+a+4)«+2|S"-i| 



Towards a contradiction assume that i3(i+a)p(0) touches flt{u) from the inside for the first time at 
<t^ <T. More precisely let 

U = inf{r > t > : S(i+„)^(o) n flt{uf ^ 0} 

where we have assumed that the set being infimized over is non-empty. Note that > by comparing 
u with a radial subsolution starting on some ball slightly larger than B(Y+a)p but still contained in fig 
and by continuity of the free boundary (Corollary 3.16 for example applies) 

B(i+a)p(0) C (u) and 3 e aB(i+,)p(0) n r*. (u). 



Now by Lemma 3.25 and 



Bp(0) CC ^t{u) for t £ [0,U+] 



we have that r2t(u) has p-reflection and therefore by Lemma 3.24 



Therefore we have 



Now, define 



S(l+a)p(0) C nt{u) C S( 

l+a+4)p (0). 
N^{N + 2)V 



A(i*) > A[B(i+a+4) 



h{x) 



|5^-i|((l + a-F4)p)^+2- 
N{N + 2)V 



2|5^-i|((l + a + 4)p)^+- 



which satisfies —Ah{x) < A(t*) where h > and h touches u{x,t) from below at (xo,<*). On its free 
boundary, due to the assumption on p and a, h satisfies 



F{\Dh\{x))^F 



N{N + 2)V{l + a)p 
|5^-i|((l + a + 4)p)^+2 



> 0. 



The barrier /i is a strict classical subsolution which touches u{x, t) from below at (a;*, t*), a contradiction 
of u being a viscosity supersolution. □ 



The following holds now due to Lemma [3. 25| and Lemma [3.26[ 
Corollary 3.28. Let u and f^o he given as above. Then VLt {u) e Sr,R for some r, R > for all t > 0. 
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4. Gradient Flow 

Now we consider the capillary droplet problem, contact line motion with volume preservation. We 
recall the PDE here for convenience, 

~Au{x,t) ^ \{t) in n{u) 

( Fv i 

ut ^ F{\Du\)\Du\ on r{u). 
We will consider the problem when 

F{\Du\) = \Duf - 1. (4.1) 

We restrict ourselves to the above choice of F because it gives the problem a simple gradient flow 
structure, although we expect similar results to hold for general F which satisfy Assumption |X] The 
viscosity solution theory developed in the previous section does not apply directly to this choice of free 
boundary velocity. Therefore we will instead deal with 

F{\Du\) = ma.x{\Du\^ - 1, M}. (4.2) 



for M > large and send Af oo to get results for (4.2 ). As mentioned in the introduction, the problem 
(P-V) with free boundary velocity (4.1 1 is formally a gradient flow of the energy, 



J{n) / \Du\^ + ([L5I 
Jn 

in the space of compact subsets of (we leave the metric unspecified for now). Above we call 

u^u[n] := argmini [ \Dv\^ :v £ H^{n), [ v^V}, 
J Jr« 

Recall that uq solves, 

( -Auix) = in n 

\ (4.3) 
[ u{x) — on dfl, 

where the Lagrange multiplier A[ri] is given by 

X[n] min{ / \Dv\^ : v £ H^{n), [ v^V}. 



Now we give the formal derivation of the gradient flow structure as motivation. One can think of the 
Cacciopoli subsets of as an infinite dimensional manifold such that the tangent space at il is 

Tn=L^{dn,\Dxn\), 



and the metric is specified as 



9nif,g) 



on 



Then suppose that fit is a smooth set valued solution of (P-V), i.e. u[rtt] is a smooth viscosity solution. 
In the geometric framework Qt is a path with velocity 



Tn,3v{t) = F{\Du\it)) = \Du\{tf^l. 
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Then we calculate 
d 
dt 



J{VLt) = [ 2Du ■ D 
= -2\{t) 



{\Du\^ + l)F{\Du\) 



2utDu ■ n + {\Du\'^ + l)F{\Du\) 



-2X{t) 



dt 



-2F{\Du\)\Du\^ + {\Du\^ + l)F{\Du\ 



(1 



Du\')F{\Du\)^-g,,Mt).v{t)). 



Ft 



Of course, viscosity solutions of (P-Vl lack the required smoothness to make the above framework 



anything but formal. Additionally, as far as we are aware, the distance associated with the metric g has 
no nice form even in the smooth case. Instead we use the pseudo-distance 



dist (fioi^^i) 



d{x, dflo)dx. 



(4.4) 



This pseudo-distance was used previously for the same problem in 17 and for the mean curvature flow 



in [h] and is motivated by the L^-based Riemannian structure formally given to the space of sets with 
finite perimeters (see 17 ). Unfortunately, dist is not a distance. For example, it does not satisfy the 



triangle inequality. We will address this deficiency in Lemma 4.2 



Our plan is to show the existence of viscosity solutions for (P-V) by constructing a gradient flow of 



the energy J^. There is a standard approach to defining curves of maximal slope for energies in metric 
spaces through a discrete time approximation. Some early works on this subject are 19 23 and the 
recent book [s] contains a quite general treatment of the method. As noted previously our problem does 
not live in a nice metric space, however the essential idea of the construction is the same. In 17 it 



is shown that the curves of maximal slope with dist satisfy a barrier property with respect to strict 



smooth sub and supersolutions of (P-Al. From the a priori estimate Lemma (3.251 any solution with 



sufficiently round initial data of (P-V) will be in Sr for all time. With this in mind, it is reasonable to 



restrict the class of admissible sets in the gradient flow to only include sets which are uniformly strongly 
star-shaped with respect to the origin. With this additional regularity it becomes possible to show that 



the solutions of the discrete gradient flow scheme will converge to a viscosity solution of ( P-V ) 



4.1. Discrete Gradient Flow Scheme. We will now define the discrete JKO type scheme to construct 



approximate solutions of the gradient flow ( P-V ) . We make a notational choice to denote these solutions 



with Lj. This is to distinguish between solutions of the gradient flow and which we use to denote 
the positive phase of a height profile. Let wq be the initial positive phase and let uq = u[uJo] be the 
corresponding initial height profile. Under the assumption that wq € Srg for some tq > 0, let us define 
a discrete scheme approximating the gradient flow of the droplet energy J" as follows: 



(4.5) 



Definition 4.1. Let M, r, h > we deflne wo^h = ujq, and define iteratively, 

i^k+i,h = arg min|j(a;) + /i^Mist (wfc^h,^) : uj e A'^{ujk,h)^ ■ 

Here A^{-) is the class of admissible sets which we will take to be 

A'\uj) {U e Sr„ and ^^(C/,^) < Mh}. (4.6) 
We suppress the dependence M, rg and h in the above notations when there will be no confusion. 



In order to show that to show that the iteration defined in (4.51 is well defined Now we apply the 
Compactness Lemma [A.4[ 



Lemma 4.1. The iteration (4.5 1 is well defined. 
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Proof. We show that for any G Sr^B, the infinium from (4.51, 

inf + /i^idSt^(wo,a;) : uj G A'^loq)} , (4.7) 

is achieved by some uj £ A'''{ujo). Let ujj E A^{ujq) be an infimizing sequence for (4.7|. In particular 
from the admissibility conditions (4.6), 

UJj G Srg.R+Mh- 

Applying Lemma A. 4 up to a subsequence the ujj converge in (Ih metric to some cu G Srg^R+Mh and we 
easily derive, 

dH{uJ,uJo) < Mh, 

I I 

so that UJ G A {uJo). Then, from Lemma A. 3 J^{-) and dist {ujq,') are continuous with respect to 



convergence in {Sro,R+Mh,dH) so w achieves the infimum in ( |4.7[ ). 

□ 

Note that since uJk is admissible at each stage we automatically get the discrete energy decay estimate, 

1 2 

J{oJk+i) + ^dist {ujk+i,0Jk) < Ji^k)- (4.8) 

We would like to use the triangle inequality (which we don't have) to get the energy decay inequality 
over time of the form: 

1 2 

— dist {uJk+n,UJk) < J{^k) - J{i^n+k)- 

nil 

It turns out that this is possible as long as we are restricting to strongly star-shaped sets. 

Lemma 4.2. Let wq, ■■■T^Jn G Sr.R, then we have the inequality, 

C(r,i?, A^)~2^ ^ ^^2^ ^ ^^^^ 
dist (cj„,a;i) < > dist (Wj+i,Wj) (4.9) 

In the following it will be useful to introduce some notation. We say that ^ < _B if there exists a 
constant C > such that A < CB. We say A ^ B or A is equivalent up to constants to B ii A < B < A. 
The dependence of the constant C on various parameters qi is expressed by writing or ~gi,....q„. 

If no dependence is expressed, then the constant is assumed to depend at most on the dimension N. 

Proof. Let uJa,^b G Sr r, then we define the distance between the sets uja and uib in the direction 

de{iOa,^b) = \Xa{e) - Xb{e)\ where {X,{e)} = duj.n {x : = 9}, i e {a,h}- 

Notice, due to the strict star-shapedness of the w,, there is only one point Xi{9) along the direction 6 in 
duji. It is easy to check that dg has the triangle inequality: 

de{uJa,UJb) < dg{uJa,^^c) + de{uJc,l^b), 

for any bounded strictly star-shaped sets Wj, i G {a, b, c}. 
Claim. The dist is equivalent up to constants to a true distance: 

disi{cJa,(^b) ^r,R,N (^J ^ ^dg{uJa,UJb)^d9^ . (4.10) 

Using the claim along with Cauchy-Schwarz and the triangle inequality for dg we will get, 
V'dist {u>j+i,ujj) >r,R.Ny^ / de{ujj+i,u>j)^dd > - / d0{ojn,ujiyd9 >r,R,N dist (a;„,a;i). 
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So it suffices to show the claimed equivalence ( |4.10 |. We calculate: 

dist (wajWfc) = / d{x,dLJa)dx= / / d{p6,dujb)luj^AiOt{p6)p^~^ dpd9 

JLOaAuJt Jo 

^r.R,N / / d{p9,dujb)'i-ij^Au,kip(^) dpd9 
■/s"-i Jo 

rmax{|X„(e)|,|X6(9)|} 

d{p6, dujb) dpdO 

IS«-^ Jmin{\X^m,\X,(e)\} 

^r.R / / \p-\Xa{e)\\ dpdO 

Js"-^ Jmm{\X^{e)\,\X,{e)\} 

p dpdO = - dg{LOa,UJbTdd. 

Jo ^ Js^-^ 

Here we have used an easy consequence of Lemma 3.10 for x S the radial distance to duji (i.e. 

— x\) is equivalent to d{x, diS) up to constants depending on r and R. In particular for x G \ 
from Figure [2] 

SO that 

T 

\x- Xi(x)\ > d{x,duji) > —\x - Xi{x)\. 

H 

□ 



As in the proof of Lemma 4.1 the restrictions on the admissible sets will allow us to take a continuum 
limit of the discrete time scheme. In particular, the assumption of uniform strong star-shapedness gives 
equicontinuity of the contact lines duj^ in space, and the Hausdorff distance movement restriction gives 
equicontinuity in time. 



Lemma 4.3. (Continuum Limit) Let uJk.h be given as in (4.5). Then the following holds for the 
interpolants 

uJh{t) = ujk,h for t e [kh, {k + l)h). (4.11) 
(i) There exists R = R(rQ,N) > such that ujh{t) C Bj^{0) for all h,t > and therefore 

Lip{dLOh{t)) < C(ro,7V). 

(ii) ujh{t) (up to a subsequence) locally uniformly converges in t ^ [0, +cx)) to uj{t), which satisfies 

dH{uj{ti),Lo{t2)) < M\ti - t2\ for all ti ><2 > 0. 



(Hi) For C — C{ro, R, N) given as in Lemma LJhit) satisfies the energy decay estimate 

^ di^t^{oJhiti),u;h{t2))<JiuJh{ti))~Jiu}h{t2))- (4.12) 



t2-ti 

Proof. We give a sketch of the proof for (i). Let x E dujh(t) and recall that 

\LOhit)\ < JiMt)) < Jiwo), 

and uJhXt) has the interior cone at x, 

{y:y-x>Oa.nd{x-y)-x>\x- y|^|x|2 - rg} U 5^(0) C uj^it) 
This cone contains ~ la^l/rg balls of radius ro/2 centered along it's axis, so that 

\x\<Riro,N) ■.= r(^''j{uJo). 



Now we prove (ii) and (iii). For simplicity we will take h ^ 2 ^ and then call the interpolant (4.11) 
ujj(t). From (i) and Lemma A. 4 the sequence ujj(t) for t G /iZ+ has a subsequence in j which converges 
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in HausdorfF distance. Diagonalizing, we get a subsequence which we continue to call j along which ujj{t) 
converges in Hausdorff distance to a limit we call From (4.6) and the triangle inequality we get: 



d/f (ajj(ii), ^^(^2)) < M\ti — ^2! for ^1,^2 G hZ'^ and for sufBciently large j, 

and uj{t) inherits the same estimate. Lo{t) is densely defined and uniformly continuous and thus has a 

unique M-Lipschitz extension to [0, +00). The local uniform convergence then directly follows. 

. . 2 

Estimate (4.121 follows from the continuity of J'{-) and dist (•,•) with respect to convergence in 

□ 



Hausdorff distance: see Lemma TA. 31 

4.2. Barrier Properties for the Gradient Flow. In |17j a discrete gradient flow is constructed for 
(P-V) without restrictions on the star-shapedness of the domain or the maximum speed of the free 



boundary. There it is proven that the discrete solutions satisfy a barrier property with respect to 



smooth strict sub and super-solutions of (P-A). An analogous result can be proven for the solutions 



of the restricted iteration scheme defined in (4.5 1. The cost of the restriction to nicer sets in (4.5) is 



that the class of barriers for which the sub and super-solution properties hold is reduced. The necessary 



conditions on admissible barriers can be deduced by inspecting the proofs in Section 3 of 17 . We give 
modified statement of the super-solution barrier property applicable to our situation below. 

Theorem 4.4. (Grunewald, Kim) (Super-solution barrier property) Let Uh — "[w/i] he the profile after 
one iteration of the scheme with time step size h > 0, and let 

A = min{A[wo], A[cj,i]}. 

Let p > 0, Xq E and suppose there exists (f) smooth with \D(f>\ ^ Q in Bp{xQ) x [0, h] with a 5 > Q 
such that 

-A(j){-,t) < X- 6 in Bp{xo)x[0,h] 



Do 



i\D<t>\^-l)<-6 



r((/.) n (Sp(xo) X [o,/i]). 



Moreover, we require that there exists C > which does not depend on h such that the sets 

LUh U ({0(-, /i) > 77} n Bpixo)), -h^<V< Ch^'^ 
are in the admissible class A]^^. Then, for < h < ho — ho{p, (j),6), the following holds: 

If^{-,0) < uo{-), then (j){-,h) < Uh{-). 



(4.13) 



Essentially, (4.13) requires that the perturbed sets of by the barrier (j) which crosses Uh from below 



stays in the admissible class. An analogous barrier property with respect to supersolutions holds as well. 

Now we show that the barrier property carries over to the continuum limit fl(t). 
Lemma 4.5. (Restricted super-solution property for the gradient flow) Let uj(t) be a continuum limit of 



the discrete gradient flow described in Lemma ^.5 Let (j) € C ' (0(0) n{0 < t < T}) such that \D(j)\ ^ 



on r(0) n {0 < t < T}, the initial strict separation fio(0) CC i^(0) holds and (j) is a strict subsolution, 
-A0(-, t) < X[uj{t)] in n{(f>) n {0 < i < T} 

(4.14) 

4>t - max{|L>(/)|2 - 1, M}\D(t)\ < on r(</>) n {0 < t < T}. 
Then (f> cannot touch u[U,{t)\ from below at any free boundary point {xQ,to) G Y{u) such that 

D<l){xo,to) ■ Xq < -ro\D(j)\{xo,to). (4.15) 



Remark 4.6. Test functions </) e C'^-^{n{(l>) n {0 < < < T}) with \D(f)\ ^ on T{(f)) can be extended to be 



C 



[o,r]). 
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Proof. Suppose not. Then there exists (j) G C'^'^(M^ x [0,r]) which touches u — u[w(t)] from below in 
M.^ X (0,r) without satisfying (4.15). From the continuity of the free boundaries of 4> and u, this occurs 
for the first time at some (xo,to) G T{u) n r(0), with < to < T. That is, there exists p > such that 

(f){x,t) <u{x,t) for e Q -Bp(a;o) X [to - p, to]- 

We may assume that (j) touches u strictly from below by making the replacement, 

(j) ^ {(j) ~ e\x - xop - e(to - t)) + . 

In particular, we will have, 

(j){x,t) <u{x,t) for {x,t) e Qnr2(0)\{(a;o,to)}, 



and for e > sufficiently small all the strict subsolution conditions for along with (4.151 still hold. 
Furthermore, we can make a subsolution which crosses u from below in Q by making a small translation 
in the normal direction to Ttg {4>) at xa , 



(p(x, t) := <f>{x + ev, t) 
iently s] 

(p(a;o,to) > = M(a;o,to) and also 



-^(2^0,^0) 



|D0| ' 



where 

Here e > is chosen sufficiently small based on the (in space) norm of (p so that (4.161 holds for <p. 



Lp{x,t) ^ u(x,t) for {x,t) ^ dpQ. 
Then we take p small depending on the modulus of continuity of Dcj) such that there exists S > with: 
\Dip\>S, tpt < [M - 6)\Dtp\, u-ip> 6 ondpQnn{ip), 



and (4.151 hold on all of Q, that is 

D(p{x,t) ■ x < -~r\Dip\{x,t) 



for {x,t) e Qnn{ip). 



(4.16) 



In order to make use of Lemma |4.4| we need to take this information back to the approximating sequence 
of discrete solutions. Let hj — be the sequence along which the discrete gradient flow scheme 

dHiLOh,{t),Ujit))^0 



uniformly for tg — p < t < tg. From Lemma A. 3 itj(-, t) := u[ujfi - (t)] converge uniformly in (a;, t) to u. Let 
j be sufficiently large [hj small) so that Lemma 4.4 will apply and also hj < d^, ip{xo,tQ) > Uj^xo^to) 



and 

ip{x,t) -< Uj{x,t) for {x,t) € dpC 
Now there exists tg — p — hj < khj < to such that 



(4.17) 



for all X £ Bp{xo) 



ip{xi,tk+i) > = Uj{xi,tk+i) for some xi € Bp{xo). 
To apply Lemma |4.4| and get a contradiction it is sufficient to show the following: 
Claim. Let 77 > —h^ then, calling Whj^k = ^^k- 

({'/'(•, ifc+i) > 77}nBp(a;o))Ua;fc+i e A^^{ujk). 



Let X £ Urp we claim that: 



d{x, ujk) < Mhj 



If a; G Wfe+i then d{x,(jJk) < Mhj from the admissibility condition (4.6) for uJk+i- If x G Qt^^iif ^ ^ 
Bp{xo) then in the case ?7 > from the containment fJt^, ((p) H Bp{xQ) C ujk and (f>t < M\D(j)\ we get, 

d{x,uJk) < d{x,ilt^{ip)) < Mhj. 
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In the case -S'^hj < -/i^ < r] < 0: 



-■q 
5 

< Mh, 



Let X ^ ujk then d(x, U.q) < d{x,LOk+i) < Mhj. Combining with the above arguments, 

dHiUr,,U!k+l) < Mhj. 



That Uri G is more subtle and this is where the condition (4.15) comes in. For a given x e dUrj, 
we need to show that for every y £ B,.g{0), 



= y + s{x~y), se [0,1]} C [/„ 



(4.18) 



This is trivially true for all x E dUfj D dcok+i from the strong star-shapedness of ojk+i- So we consider 

X e dUr, n d{nt,+r i^-v)r\ Bp{xo)) c Bp{xo) n ar^t,^, i^-v)- 

Let y € (0) and we consider the line segment from x to y parametrized by 

z{s) = x + s{y — x). 

First we show that - for s > small depending on the norm of - z(s) e C/^. Using \y\ < r and 
( p6t , 

0(z(s)) = + sDv3(x) • (y - a;) + o(s) > s(-|y| + ro)|D<^|(a::) + o(s) > 0. 
Suppose towards a contradiction that z(s) exits C/^ before s = 1, then by continuity it must pass through 
dUr, for the first time at < sq < 1 defined by, 

1 > So = inf{s > : z{s) G dUr,} > 0. 
If z{so) e dujk+i then we are done due to uJk+i G 5^: 

Lz{so),v ^ C [/,,. 

So it remains to consider the case z{sq) € -Bp(xo) H 951fj.^j ((/?) and in this case: 



0>-^(z(.)) 



Dip{z{so)) -(y-x) ^ {I- so)Dip{z{so)) ■ (y - z{so)), 



rearranging this and using |y| < tq we get 

D(p{z{so)) ■ z{so) > D(p{z{so)) ■ y > -ro\D(p{z{so))\ 



which contradicts (4.16). 



□ 



An analogous statement about touching from above by strict classical super-solutions holds as well. 
Lemma 4.5 almost says that uj{t) is a viscosity super-solution of (P-V) except for the restriction (4.15) 
on the barriers. In fact we can make the following statement: 

Proposition 4.7. // there exists r' > ro such that uj{t) G Sr' for all t e [0,T] then u[uj{t)] cannot 
be touched from below (above) by any strict classical sub-solution (super-solution) of (P-X[uj(t)]) and is 
therefore a viscosity solution of (P-V). 



Proof. This is simply because any <j> £ C^'^(E^ x [0,T]) which touches from below (above) at a 

free boundary point must necessarily satisfy (4.15). To see this fact, let x £ duj{t), now from the strong 
star-shapedness of i^(i), 

{x + s{x-y) ■.s>0 Sindy£ Br'{0)} C R^\uj{t) C ]R^\17t(0) 

and so for any y £ B^.' (0) 

d 



0> —(j){x + s{x-y)) 
as 



D(j){x) ■ (x-y). 
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Rearranging above inequality, we obtain ( 4.15[ ) 



D(l){x)-x< inf D(l){x)-y = -r'\D(l)\{x)<-ro\D(l)\{x). 



□ 



Unfortunately this statement is not that useful to us because we would still need to show that the 



gradient flow preserves strong star-shapedness. That is if we had a Lemma for uj[t) like Lemma 3.20 
then from the above Proposition we would know that uj{t) was a viscosity solution of (P-V) at least 



for a short time, and then we could apply a result about preservation of strong star-shapedness for 



viscosity solutions like Lemma [3.26 to show that uj{t) is a viscosity solution globally in time. Instead of 



attempting to prove Lemma 3.20 for solutions of the gradient flow we will show directly that uj{t) and 



the viscosity solution of (P-A[a;(t)]) are the same. The idea comes from the following corollary of the 
proof of Proposition |4.7[ As usual the corresponding result for supersolutions is also true. 



Corollary 4.8. (Comparison with star-shaped classical suhsolutions) Suppose (p e C^'^(] 



[Q,T]) IS 



a strict classical sub-solution of (P-X[uj{t)]) and the initial strict separation 

holds. Moreover for some r' > ro the positive phase of (j), U,t{(j)), is in Sr' for all t G [0,r]. Then 

0(-,r) <u[c^(T)]. 

4.3. Comparison. Let Lo{t) be the continuum limit of the discrete gradient flow scheme defined above. 
We want to show that u = 7i[w(t)] is a viscosity solution of the problem, 



-Aw(a;,<) = \[uj{t)\ 



u{t) 



= min{|L>i;|2 - 1, Af}|Du| on du{t), 



(4.19) 



at least as long as the unique viscosity solution of (P-A[w(t)]) satisfies the appropriate strong star-shaped 
condition. 



In fact u satisfies a comparison principle with respect to sub and super-solutions of (4.191 which are 



star-shaped with respect to a slightly larger ball than Brg (0) . This is the non-smooth version of the idea 
of Corollary 4.8 Various iterations of the proof can be found in the literature [6l[T0 20 , we will only 



sketch most of the details of the proof the essential idea is, as in Corollary |4.8[ to demonstrate that the 
restricted barrier property of the gradient fiow solution Lu(t) is sufficient to make the argument work. 



Proposition 4.9. (Comparison with sub-solutions) Letu{-,t) 
and V a continuous sub-solution of 

-Av{x,t) < X[uj{t)] 



u[Lu(t)] as above from the gradient flow 



with initial data ordered. 



Vt < min{|L»u|2 - 1, M}\Dv\ on 



na{v) C io{0). 



Let r > from the restriction on admissible sets in the discrete scheme (4.5 1. Suppose there exists 
r' > To such that ilt(w) G 5,' for all t G [0,T]. Then 

v{-,t) < u{-,t) on [0,T]. 

Proof. 1. First we show it suffices to treat the case when u and v are initially strictly separated, 

na{v) CC uj{0), (4.20) 

and there exists e > so that v satisfies, 

-Av{x,t) <{l-e)X{t) for {x,t)en{v). (4.21) 
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Let us assume temporarily that result holds in this case and we show that it also holds under the non- 
strict conditions given in the statement of the proposition. This is accomplished by a similar device 
to the one used in the proof of the strong comparison type results Lemma |3.12| and Proposition |3.11| 
Let e > small enough that (1 + e)~^r' > tq, c from Assumption [b| and < a < ro then define the 
perturbation, 

=sup|,|<,,_,,t(l + e)-3w((l + e)(.T + z),<) for < t < f . 



Now the perturbed subsolution satisfies (4.20) and (4.211 and so we get, 

v{x,t) = swp^yQV^{x,t) < u{x,t) for 0<t<a/c 
Iterating this [Tc/a\ times we get the desired result. 



2. Now we work in the case when (4.20) and (4.21) hold. Suppose towards a contradiction that v 
crosses u from below on X [0,T]. In order to have some regularity at the touching point we use the 
space-time convolutions described in Appendix [b| Let e,6 > with 6 < e/T define the sup-convolution 
of V 

V{x,t)^ sup v{y,s) (4.22) 

\iv,s)-{x,t)\<e 

and the inf-convolution of u 

U{x,t)^ inf Z{y,s), Z{y,s):= \ni u{w,t). (4.23) 

\{y,s)-(x,t)\<e~St \{-w,T)-(y,s)\<e 

defined in the domain 

:=M^ X [e,r-e]. 
We will show that V still has the following properties: 

(a) Vlt{V) strongly star-shaped with respect to a ball of radius larger than r, 

(b) y is a subsolution of (P-A[a;(t)]) and is initially strictly separated from [/, 

^,{v) cc n^{u). 

First we show (a), let e < (r' — r)/2, then 

^t{V)= U y + nt+s{v). (4.24) 

\{v,s)\<e 

Noting that all of the sets in the above union are star-shaped with respect to B(r+r')/2(0) implies that 
^t{V) is as well. This proves (a). 

Now we show (b). Choosing e smaller if necessary based on 77 and the modulus of continuity of \{t) 
the subsolution property for V is from Lemma |B.1[ Then simply from the initial strict separation of u 
and V along with the Hausdorff distance continuity of the free boundaries of u and wwe also have the 
strict separation of U and V at time e 

a(v") cc n^{u). 

Now by our assumption we know U crosses V from below for the first time at some point Pq = 
(ccoj^o) G ^e- By the strong maximum principle for subharmonic functions 

xo^Tt,{U)C^Tt,{V). 

At the touching point Pq the positivity set of V has an interior space-time ball of radius e centered at 
some point Pi = (xi,ti) G r(u), i.e. |Pi — Po| = e and 

i?i := {(y,,s) e X [0,r] : |(y,,s) - Pi| < e} C n+{V). 

This ball has the tangent hyperplane through Pq 

i7i = {(y, ,s) e X [0, T] : ((y, s) - Po) • (Pi - Pq) = 0}. 
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Let (yTfri) be normal to Hi in the direction inward to ^{V) with \v\ = 1, m is the advance speed of 
the free boundary T{V) at Pq- We will show below that m is finite, but at this point let us include the 
possibilities that m = oo orm = — oo. 

Similarly, there exists P2 = (3:2,^2) G r(Z) with IP2 — Pi\ = e — ^to and fl{U) has an exterior 
space-time ellipse at Pq of the form 

B2 := {{y, s) G X [0,r] : |(y, s) - F2I < e ~ Ss} 

while there exists (^3,^2) G r(u) with \x2 — x^l = e and Q{u) contains the following set: 

E = {{y,.s) eR"" X [0,T] : \y ~ w\ < e,\iw, s) ^ P2\ < e - 6to} d n{u), (xg, ^2) € ^ n r(u). (4.25) 

Let H2 be the tangent hyperplane to B2 at Pq and i/2 be the tangent hyperplane to B2 at P2- 

Lemma 4.10. is not horizontal. 



Proof. We refer to the proof of Lemma 2.5 in 20 , which rules out the possibility that m = +00 



Next let us rule out the possibility that m — —00. Let X3 and E be the set given in (4.25). Due to 
the star-shapedness and Holder continuity of Q{u), it follows that {x3,t2) lies in the spatial boundary of 
E n{t ~ t2}. We consider the classical subsolution 

Take t > sufficiently small that Dq :— E D {t — t2 — t} is nonempty. Let Di be the top portion of 
the closure of E, i.e., Di :— E O {t = t2}, and let us define the interpolation of Dq and Di, i.e., 

Dt (1 — s(t))_Do + s{t)Di, where s{t) is linear and satisfies s{t2 — r) = 0, s{t2) = 1- 

We choose r small that s'{t) < —1. Now let us consider the space-time domain 

S :^ U m) ^D{t)] 

T<t-t2<0 

Then, due to the fact that T, C E, the positive set ^l{u) crosses S for the first time at (^3, t2). Let us 
consider the classical subsolution tp in the domain satisfying 

-A^(f){x,t) = X[w(t)] for {x,t) e E 

= for {x,t)edDit) ^^26) 



= mm 



> for {x,t) e ^dD{t) 



Then ^ crosses u from below at (2:3, ^2)- One can check from the definition of u and the strong star- 
shapedness of n{V) that the outer normal v of Di at X3 satisfies u ■ X2 < — tq and thus it contradicts 
Lemma 14.51 □ 



Let (v' , m') be the normal to the hyperplane H2 rescaled as before. Then as in Lemma B.l the normal 
to H2 is {v' ,m' + 6). As a consequence of the ordering V{x,t) < U{x,t) for all t < we also get the 
ordering of the advance speeds, 

m' + 6 < ni. 

Since U <V for t <tQ we get the following inclusions 

Si n {< < to} c {t/ > 0} n {v^ > 0} and i?2 n {t < to} c {J7 = 0} n {V = o}. 

In particular ^lt„{U) and VLtoiV) both have interior and exterior spatial balls at xq, 

Bir\{t = to} and i?2 n {t to} 

which must both have centers lying along the same axis, or equivalently both free boundaries have spatial 
inward normal v at xq. From the strong star-shapedness of fitjj(y), and this is a key point as we have 
seen in the proof of above lemma: we get that 

vxo< -r" < -ro. (4.27) 
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Now the idea is to construct a smooth strict sub-solution (f) which touches y{V) from below at (a;o, io)- 
Then a translation of will touch u from below at P2 from the inadmissible direction, 

(2^0,^0) 



leading to a contradiction of Lemma 4.5 
Lemma 4.11. 



-1 < m < M. 



Proof. The fact that m < M follows from a relatively simple barrier argument, based on the condition 
Vt < MiZJwl on Tt{v) as well as the fact that ri(w) has an exterior space-time ball at {xi,ti). 

It remains to show that m < —I. If it is, then 

m' < -1 - 6. (4.28) 

Take t > sufhciently small that Dq := B2 d {t = t2 — t} is nonempty. We can now construct a classical 
subsolution ip in the domain 

E := Ut2~T<t<so{B2 - 0^2), 



similarly as in the proof of m > —00 in the lemma above, and use (4.28) to derive a contradiction. □ 



Lemma 4.12. Near the point Pq we have the nontangential estimate 

. V{xo-dv,tf^) — 

™ ; > \/in + 1. 



lim inf ■ 



(4.29) 



Proof. The proof is based on construction of the barrier for v to yield a contradiction, in the event that 
the lemma holds false, and it is parallel to the proof of Lemma 2.6 in ,20j . 

□ 



Now as in 20 for any 77 > we construct a smooth test function 4>{x, t) 
following properties: 



I a; — XQ\,t) with the 



-A^<j>{x,t) < 


for 


{x,t) ^ \B2 


> 


for 


{x,t) e B2 


< 


for 


{x,t) ^ B2 


\D(l)\{x,t) = 


- 1(1 — r;) on 


dB2 n {t = ta} 



(4.30) 



Now from the definition of m' we have for x G dB2 n {t = ta}, 

T^^{xM) = mf <m-5 

and choosing t sufhciently small depending on the norm of (j) and rj small depending on S and m we 

get for x e dB2 n {^3 - t < t < ta}, 

^(x,t2)<m-<5/2<|D,^p-l. 
Therefore (/> is a strict subsolution in the region 

{B2 \ 3B2) n{t2-T<t< t2}. 

Next we show that touches u from below at P2 — {x2,t2). For d > sufficiently small depending on 
the norm of (p we have 

(j>{x, t) < VmTT(l - 2ri)d on d{l - d)B2 n {ia - r < t < ia}- (4.31) 

Let (x, s) e B2 and let 

d := d{x,dB2 n {t = s}) 
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then since Pq is the center of B2 wc have that |Po ~ {x, s)\ = e — 5to — d and 

|(xo + dv,to) - {x, s)\ < \Po - {x, s)\ + d ^ e - Sto. 



Therefore from the definition of U as an infimum and from Lemma 4.12 for any 77 > there exists do 
such that d < do imphes 

Vm + l{l - 2r])d < U{xa + dv, to) < u{x, s). (4.32) 
Now combining ( 4.31[ ) and (4.32) with the fact that B2 C and (j) — Q on dB2 we get that 

(l){x,t) <u{x,t) for {x,t)e{d{l-d)B2yJdB2)C^{t2-T <t<t2}. (4.33) 
Then since u — (p is superharmonic by the strong minimum principle we get, 

(t){x,t) <u{x,t) for (a;,t) e S2 \ (1 -d)B2 n {^2 -T < t < ^2}, (4.34) 
= <t>(P2) = SO (f> touches u from below at P2- This is a contradiction of Lemma |4.5 



and ^(^2) = 
from (|4.27|) 



D<j> 



{X2, 12) ■ X2 = V ■ X2 < V ■ Xo < -rg. 



Since 



□ 



5. Convergence for Solutions with Global in Time Star-Shapedness 

Now we can combine the results of the previous sections to get our main result. Under the assumption 
of sufficient roundness of the initial data JIq phrased in terms of p-reflection any continuum limit of the 
restricted gradient flow of the functional J is also a global in time viscosity solution of the problem 



(P-V). To make this precise let us define the class of weak solutions arising from the gradient flow 



scheme described in section |4l 



Definition 5.1. An evolution uj{t) : [0,oo) — 5,- is an energy solution of (P-V) if there exist aJA4(t), 
Mk — > 00 as /c — >■ 00, which are minimizing movements of the restricted gradient flow scheme from 
with initial data a;(0) such that, 



Definition 



4.1 



A droplet profile u 
wit). 



d/f (cjMfc (^), ^^(^)) ^0 as A; — > 00 locally uniformly in [0,cxd). 

: X [0, 00) — > [0, 00) is an energy solution if u = 7i[w(t)] for an energy solution 



Now we can say the following about any energy solution arising from an initial data with p-reflection. 
Theorem 5.1. Let V > and ilo a domain in such that ilo has p-reflection with p satisfying, 

p < CnV"^ , 

where is a dimensional constant. Then there exists an energy solution u ; x [0, -fcxi) — > [0, +00). 
Moreover any energy solution u, the following holds: 

(a) u is also a viscosity solution of the free boundary problem, 

—Au{x,t) — X(t) in rttiu) 



ut = {\Du\'^ - l)\Du\ on Ttiu) 

^ {u(-,0)>0} = f^o, 
where X{t) is chosen to enforce the volume constraint for all t > 0, 

u{x, t)dx = V. 



(P-V) 
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(b) The positivity set flt{u) has p-reflection for all t > with iT^ixert{u)i\x\ ~ p) bounded below. The 
energy J is decreasing along the flow and additionally u satisfies the energy decay estimate for 
allt> s> 0, 

^d5t(r!,(u), rit(u)) <p,^v J^^siu)) - J{nt{u)). 

(c) The flow of the sets fit(u) converges uniformly modulo translation to the radially symmetric 
equilibrium solution, 



2N 



where r^ is given in (2.2 I and A* can be calculated from the volume constraint. More precisely 
we mean that, 



mfldni^tiu), Br*{x)) : x G Sp(0)} as i oo. 

Proof. 1. The existence proof as well as the energy estimate follow from a fairly straightforward com- 
bination of the results we have already proven. We give an outline of the proof. Given the assumption 



on f2o having p-reflection we know due to the apriori estimate Lemma 3.26 that any solution of (P-Vl 



will be in Sr for some r > as long as it exists. Then we construct the restricted gradient flow solution 
with initial data fig as in Section [4] where the restriction on the radius of the strong star-shapedness is 
strictly weaker than that from the apriori estimate. Given the Lagrange multplier \{t) associated with 



the restricted gradient flow we then solve (P-A). Then the idea is that the viscosity solution of (P-A) 



is strongly star-shaped with a larger radius than the restriction on the gradient flow and so we will be 



able to use Proposition 4.9 to show that the two solutions agree for all time. 



Let 6 > so that 



5(i+b)p(0) C n and 



(1 + b)p < CnV«+^ . 

b)p will be contained in ^It for all 



Then we expect thanks to the apriori estimate Lemma 3.26 that (1 
time. In particular we expect Vtt will be in Sr for r = bp. Let R — R{r) > be so large that any set in 
Sr which touches Bj^ from the inside must have larger energy than VLq and then define 7 = ^ • That it 
is possible to choose R in this way is described in Lemma 4.3 Let A/ > and let uJi\i{t) — uj]\i^jr{t) be 



the restricted gradient flow solution with initial data fip as constructed in the last section. The notation 
indicates that LOM,'yrit) is restricted to remain in Sr and its free boundary can move no faster than M. 
We suppress the dependence on jr. Define the putative Lagrange multiplier Xm to be 

AmW := AKf(i)], 

and let um be the possibly discontinuous viscosity solution of (P — Aj\/) constructed by Perron's method 
in Theorem 3.15 We will show that um and are the same. 



Let / be the largest interval containing the origin on which um and zi[a;A/] agree. We will show that 
/ is open and closed in [0, -|-cxd). Since um agrees with on / it is continuous and has constant 

volume and thus it is a volume preserving viscosity solution. Therefore Lemma |3.26| implies that: 

^^i('"Af) has /9-refiection and -B(i+6)p C Qtiunj) for t G /. 

In particular f2((uM) is in Sr on /. Suppose / — [0,T) for some T > 0, then the short time existence 
theorem implies that, for some small tQ{r,R), um is continuous on [0,T -I- to). Since the set where two 
continuous functions agree is closed / = [0,T]. Suppose / = [0,T] where T may be equal to zero. Then 
thanks to Lemma [3.20| there exists tM > such that: 

^t{uM) e Sr' for some r' > jr on [0, T + t^f). 



Then from Proposition 4.9 m[wm] satisfies a strong comparison principle with respect to viscosity so- 
lutions which are strongly star-shaped with a larger radius than jr so um = u[^m] on [0,r -I- tj^j). 
Therefore / = [0, +oo) and there exists a global in time continuous viscosity solution of (P-VUf which 
has /9-reflection and is in Sr,R. for alH > 0. 
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2. Now we show the existence for (P-Vl without the restriction on a maximum speed. The key point 
in this case is the equicontinuity afforded by the fact that the ilt(ujv/) G Sr^R with r,R independent of 
M. In particular from Corollary 3.16 we get for some C > and a e (0, 1) independent of M, 



dH{rt{uM),rs{uM)) < c\t- 



(5.1) 



Then we also derive thanks to Lemma A. 3 the equicontinuity of the Lagrange multipliers AAf(i)- Taking 
a subsequence such that the AA/(i) converge uniformly on compact subsets of [0, +oo) to some X{t) we can 
apply Lemma |XT9] We get that along this subsequence the um converge uniformly to a viscosity solution 
u of (P-Vl with free boundary velocity F{\Du\) = \Du\'^ — 1. Due to the uniform convergence of ftt{u]\f) 
to f2f (um) in hausdorff distance sense the energy estimates for the um and the Holder regularity in time 
(5.1 ) carry over to u. Then u is an energy solution by the definition. Unfortunately by the compactness 
method we do not know whether there is any uniqueness of the limiting Lagrange multiplier X{t). 

3. Now we show that any subsequential limit of the u{-,tn) must be a viscosity solution of the 
equilibrium problem (EQ|. Note that the same result is true for all the um- 

Claim. Let — > oo such that Q,t^{u) converges in Hausdorff topology to Q,oo- Then u[rioo] is a 
stationary solution of (P-V), that is it solves the equilibrium problem (EQ) in the viscosity sense. 



Define [/„ : [0, +cxj) — > S^^r by 



Un{t) ■■= n 



t+t„ 



and we also consider the viscosity solutions of (P-V) which lie above the Un(t), 

Vn{x,t) :— u[Un{t)]{x) with Lagrange Multiphers A„(i) := X[Un{t)]. 

First we show that, uniformly in i > 0, 
Since is J{Un{t)) is monotone decreasing for all n > we have for alH > 

J{n^) - inf J{^s) < J{Un{t)) < J{UrM) = Ji^tJ, 

but due to Lemma [A.3| along with the convergence r2„ — > Oqo in Hausdorff distance, 



(5.2) 
(5.3) 



Now we will show that up to a subseqeunce the w„ converge uniformly on compact time intervals. 
Recalling the uniform Holder estimates from (5.1 1 for flt{uM) which carry over in the limit to fit(u), 

dj/(C/„(i),C/„(5))<qt-5|", (5.4) 

and from Lemma [4?3| the energy estimates, 

-^dirt(C/„(s), Ur.it)) < J{Ur,{s)) ~ JiUnit)). (5.5) 

t — s 

The paths J7„ are a sequence of equicontinuous maps into (5r,_R, dh) and so from the Compactness Lemma 



A.4 



there exists Uoo '■ [0, +00) 

U„ - 



->■ Sr.R with Uoo{0) = r^oo such that up to taking a subsequence, 
Uoo uniformly on compact subsets of [0, +00). 



Now we show that Voo is a stationary viscosity solution of (P-V). From Lemma A. 3 we get the 
following: 

(i) Vn — > Voo uniformly in (x, t) on compact time intervals and therefore from the stability of viscosity 
solutions under uniform convergence - Lemma 3.6 - w is a viscosity solution of (P-V). 

(ii) dist(C/„(s), Unit)) — > dist([/oo(s), C/oo(^)) uniformly on compact subsets of [0, +00) x [0, +00). 
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Combining (ii) with (5.5) we derive the energy estimate for C/oo(') for all t > s > 0: 

-^ddst{U^{s),U^it)) < J{U^[s)) - J{U^{t)) = 0. (5.6) 
t — s 



So Voo = w[r2oojis a stationary viscosity solution of (P-Vl. Then due to Theorem 2.1 Theorem 2.2 
and Corollary 2.3 it follows that i7oo = -Br*(z^o) for some point xq g M^. Actually xq is not completely 
arbitrary since we know that fioo must have p-reflection. One can easily check that this implies € 

4. Finally we show that the convergence is uniform modulo translation. Suppose that there exists a 
sequence of times t„ — >■ oo and a 5 > such that 

i nf dH{^t„i.u),Br*{x)) > 6. 
X e-Bp(o) 

By taking a subsequence of the tn we may assume that Q,t^ [u) converges in Hausdorff distance to some 
Oqo- By part 3 of the argument u[0oo] must be equal to Br*{xo) for some xq € Bp{0). Choosing n 
sufficiently large so that 

d{nt^{u),B„{xo)) < S 

we derive a contradiction. 

□ 

Note that above theorem leaves the possibility that the drop oscillates between a family of round 
drops with its speed going to zero but not fast enough to converge to a single profile. Below we show 
that, if the drops are sufficiently regular at large times, then this does not happen. Such regularity, when 
the time is sufficiently large so that the profile of u is sufficiently close to a round one is suspected to be 
true in the light of existing results introduced by Caffarelli et. al. (see e.g. the book [o]), but verifying 
this for our setting would be highly nontrivial and thus we do not pursue this question here. 

Proposition 5.2. (Conditional uniqueness of the limit) Suppose additionally that u{-,t) are uniformly 
C^'°', then ilt — Br*{xo) for some xq. 

Proof. Let i„ — > oo be a sequence of times along which Tt^ (u) converges in Hausdorff distance. The limit 
is a ball dBr*{xo) for some xq which is compatible with the reflection symmetry of fl. By translating 
we may assume that xq = 0. Note under this translation flt{u) no longer need have p-reflection or be 
strongly star shaped with respect to the origin, this will not affect the proof. From Lemma [A. 3| we get 
that u{-,tn) converge uniformly to 

Due to the assumption, we have 

sup||M(-,i„)||ci.° < Ci < +00. 

n 

In particular the I?u(-,t„) are uniformly bounded and equicontinuous so they must converge uniformly 
to Dueq- Let (5 > small enough that Ci(5" < 1/2 and choose n sufficiently large that 

dH{Tt„{u),dBr*) < 6 and \\Du{-,tn) - DueqWoo < S. 

Now let X £ Tt^iu) and let y — y{x) £ dBr, such that \x — y\ ^ d{x,dBr,). We calculate, 

\Du{x,t„)-DuEQiy)\ < \\Du{-,t,,) - DueqWoo + Cidix,dBr,r < CS" 

so since \Dueq\ = 1 on 95^* we have that \Du{x,tn)\ > 1/2 and 

|^(a:,i„)~^|<C<5". 

Rephrasing the above in terms of the interior normal field i>n to Fj^ (u) we get that 

{i^nix),x)>l-o{l). (5.7) 
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In particular this means, thanks to Lemma |A.H that for any p > there exists n sufBciently large so 
that flf^ (u) has p-reflection. Since p-reflection is preserved under the flow for p small this means that 
for every p > there is a T such that for t >T, flt{u) has p-reflection. Thus any subsequential limit of 
the ftt(u) has p-reflection for every p > and must be a ball centered at the origin. 

□ 

There is one nontrivial case where we can say that the limiting equilibrium solution is unique. When 
the initial data is symmetric with respect to N orthogonal hyperplanes through the origin in addition 
to all the conditions given in the statement of Theorem |5.1| then the limit is unique. In this case the 
reflection symmetries are preserved by the equation and so any Hausdorff distance limit of the flt{u) 
will share these symmetries. Then it is basic to check that the only ball of radius r* which is symmetric 
with respect to N orthogonal axes through the origin is in fact centered at the origin. We record this 
fact in the following corollary: 



Corollary 5.3. If fla .satisfies all the conditions of Theorem 5.1 and additionally is symmetric with 
respect to N mutually orthogonal hyperplanes through the origin then any solution u of (P-Vl with 
initial data Q,q constructed from the discrete gradient flow scheme of Section ^ satisfies, 

dH{^t{u),Br*{Q)) as t ^ oo. 

Appendix A. Geometric Properties 

Let n be an open bounded domain in which is strictly star-shaped with respect to 0. Let p > 
such that Bp{0) C il and let H he a. hyperplane in such that H n -Bp(O) = 0. Let be the open 
half-space of H which contains -Bp(O) and be the interior of its complement. Then define 

n+ = nnH+ and = r2 n 

and define the reflection through H by 

^ X - 2{x - y, vh)vh 

where i>h is the unit normal to H (pointing inward towards i7„ for concreteness) and y Cz H . 



Lemma 3.23 , Suppose ft has p-reflection, then 



Proof. Take xq £ dfl such that 



sup \x\ — inf < 4p. 



X^OU 



Let H be any hyperplane tangent to -Bp(O) such that xq G H+. 

1. Claim: The reflections 0//(xo) with respect to the hyperplanes described above cover all directions 
in the sphere, more precisely 



^'"[^"l, : H tangent to B„{0) and xq e H+} = S^'K 
\(Ph{xo)\ 



We index H by its normal vector, so for G ^ let H^j be the hyperplane through normal to and 
then define: 

Let Lu £ S^^^. Without loss, by changing coordinate names, we may assume that xq — e\ and u) G 
span{ei,e2}. We restrict ourselves to hyperplanes with normal direction v G span{ei,e2} and thereby 
reduce to the case N — 2. Let ip such that cos ip — p/|xo|. Then for any > ip the plane through 

p((cos6')ei + (sin 0)62) 

has Xq in the same half-space as Bp. Let us consider the continuous map / : [0, tt] — > [0, tt] defined by, 
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We show that /(O) — and /(tt) — n and thus / maps [0,7r] onto itself. We chose (f so that xq G 

{xn - pj/(0),i^(0)) = {xo,v{0)) - p = |a;o| cost/? - p = 0. 
Therefore (t)u(a){xQ) — xq and /(O) = 0. Meanwhile since I'^n) — t: we compute directly that 

0iy(7r)(a;o) = (/O - ko|)ei 
and thus /(tt) = tt. Now in order to show that uj is hit we just choose 9 so that: 

=cos-i(a;,ei) 

or in other words, 

2. Let u! S S^~^, by $7 strictly star-shaped with respect to there exists t > such that toj e dil 
and s < t implies sui G fl, s > t implies sui G fl'^ . We want to show that 



t< |xo|+4p. 

I0h(2:o)| 



Let H be such that xq G i7+ and ."^"[^"j, = uj, which is possible by part 1 of the proof above. Then by 



(J 

star-shapedness Aa;o G for A G (l,+oo) and the analagous statement for the reflected domain, 

(J 

X(I)h{xq) e (I)h{^+) for a G (1, +00). 

Moreover since il has p-reflection fi_ C 0/^(17+) so in particular 

t 



7(j)H{xo) =tuj G C 4iHi^ + ) 



\4'Hixo)\ 

SO t < \cj)H{xo)\- Now, noting that —pvH G H due to H being tangent to -Bp(O), 

<Ph{xo) = X - 2{xo + pvh,vh)^'h = (a;o - (xq, vh)i^h) - (2r + {x(3,i'h))vh 
so by Pythagoras 

\4'h{.xo)\^ = |a;o - (xo, + (2r + (xo, j^h))^ 

= |a;o - (xo, -f- 1(2:0, i^_H->P +4p^ -I- 4p(2;o, i^/f) 

= |a;op-t-4p2-f4p(xo,i^H) 
rearranging (and noticing that (xq, i^ij) > — p due to our assumption that x^ G 

p2 + p(xo, I/ff) 



i < \<I>h{xo)\ < \xo\+4. 

< \X0\+4: 



\<I)h{xo)\ + \xo\ 
p'^+p\xo\ 



\(I)h{xo)\ + \xo\ 
2r 

- I^ol + 4 — 

1 + \(l)H(Xo)\/\xo\ 

< |a;o| +4p, 

where we have used |xo| > p and |(/)h(2:o)| ^ l^^ol in the If^st two lines. This completes the proof. □ 



Lemma |3.24[ Suppose has p-reflection. Then fl satisfies the following: 

(a) for all x G dQ there is an exterior cone to fl at x, 

X C(x, 0a;) C \ f7 where cos(/)^ = , G (0, 7r/2), (A.l) 

and C{x,4>x) is the cone in direction x of opening angle as defined in (3.8); 

(b) il G 5r where 

r = r{p, inf N) = ( inf |x|2-p2)i/2. 
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Proof. Let xq G ^^l and let TZ be the collection of planes which pass through xq and are admissible for 
reflection i.e. 

n = {H hyperplane in : xq € H and H n Bp{0) = 0}. 

Notice that if G 5 then for any a > the plane H + avn C iJ_ so H + avn is also admissible for 
reflection. We can also characterize TZ as hyperplanes H C such that infygjy \y\ > p and xq G iJ, 

p < inf |y| = |(2;o,i^H)|- 

So if we think of these planes as being indexed by their normal vectors, 

TZ — {H : {vh, Xq) > p and xq G H}. 
Now let J/ G + C[xQ,(f)xo) as above, then from the definition of (f)^^, 



'\y - xa\ 



,Xo) > p 



so the plane Hy with normal i^Hy = \y-x"o\ through xq is in 7?.. Now we will use the reflection comparison 
with respect to the plane 



H ^Hy + -{y-xo) 



1 

which is admissible for reflection by the remark above that H C Hy^. Then from the p-reflection 
property of il, 

y = (t)fj{xo) G \ n H+) cR^\{nr\ h.) 

and therefore y G \ f^. Since y was arbitrary, 



This proves (a), (b) follows from (a) and Lemma 3.10 □ 



Conversely we can show that given a domain 51 which is uniformly close to a ball around zero and 
has some uniform condition on the directions of its normal vectors (i.e. star-shapedness with respect to 
a large ball) has p-rcflection. So p-reflcction is a natural condition to describe closeness to being round. 

Lemma A.l. Let fl G So, n{x) be normal to dQ, and let p > such that Bp{0) C fi. Suppose that for 
H^^^ almost every x G dVl 

|(n(x),x)p>|xp-pV5 (A.2) 

and also, 

sup \x\'~-\y\' <p^ (A.3) 

then U, has p-reflection. 



Remark A.2. The conditions (A.2 1 and (A.3 1 should be interpreted together as a smallness requirement 
on the Lipschitz norm distance between dQ and the nearest sphere or alternatively that 57 G Sr,R for r 
sufhciently large depending on R and p. 

Proof. We first prove the result when dfl is C^. Then derive the result for nonsmooth fl by density. Let 
n{x) G C(M^;M^) be a vector field normal to dn. Let v G 5"^"^ and define 

H ^ {x : X ■ ^ 0}, H+ = {x : X ■ ly > 0}, ^ {x : x ■ < 0} 

the hyperplane normal to through the origin and corresponding half spaces. For s > define the 
translates of H, 

H{s) = H + sv, H+{s) =H+ + SI', (s) = + sv. 
We want to show that for s > p 

^His)mnH^{s)cnnH^{s). (A.4) 
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We make the following notations for simplicity, 

ris := fin H^{s), and fls := (t)H{s){^) H^{s). 



From ri C Bfj{0) wc know that H{R) does not intersect fl and therefore (A.4) holds for s > R. Now 
move the plane H inward towards the origin until it touchs f2 for the first time at 

So inf{.s > : H{s) n = 0}. 



Then the containment (A.4) holds until Smin('^) which we call t for convenience. We want to show that 
t < p. Note that t < sq since for s > sq the intersection ilr\H-{s) is empty and so (A.4) holds trivially. 
Moreover t < so because dV, is C^. At t there are two possibilities. The first is that fit touches fit from 
the inside at some point off of H{t), that is there exists x £ dfl D H+{t) such that: 

<t>H{t){x) e ^ntr^^'^t\H{t). 

The second is that H{t) intersects 951 at a point where v is tangent to dVl. 

In the first case, call y = 4>H{t){x) to be the point in dflt D dQ,t \ H{t). Note that (j)H{n{x)) is the 
normal to Oil at y. Initially we suppose that 

Then a simple calculation, or some geometry, shows that: 

Now we derive a bound for t from above in this situation, 

= \x- + \y\^ + 2{x -y,y)^{x~ y, + \y\^ + 2\x - y|(i/, y) 
= (2t - 2(y, v)f + |yp - 4(i - (y, v)){v, y) 
= At^-m{y,y)+^u,y)^ + \y\'' 

>t' + \y?. 

Keeping this in mind we now work in the case when 

{y,v) > I- 

Because x lies to the in H^[v) and we have assumed a bound on the angle between x and n{x) we also 
get a bound on the angle between n{x) and v, 



{n{x),v) = {n{x) -x,v) + ^>^2\l- Jl- + ^> * ^ 



Meanwhile, 4>H{n{x)) inherits the opposite bound, 

t p"^ 

{(f>Hin{x)),iy) = -{n{x), v) < -— + 
Due to our assumption on the angle between y and v we also get a bound in the other direction, 



bH{n{x)),i^) = {(l)H{n{x)) -y,iy) + (y, i') > 



5|y|2 4|y| 

Combining the above two bounds we get that, 

p'^\x\ p2 

then calling 7 = and using p < min{|a;|, \y\} we rearrange to get, 

+ ^ ^ f P ^ , P \ ^ 
t<T\ PO i I 1 +Pri - P- 
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In the second case, let x E dil D H{t) such that, 

{n{x),v)={), and therefore \{x,v)\^ + \{x,n{x))\^ <\x\^ , 
Rearranging and noting that x G H{t) impHes \{x, v) \ = t. 



\x\' - p'))Y" = p. 



Finally in the case of general Q. E Sq one can approximate dO. by boundaries of domains which 
converge in Lipschitz norm (in the appropriately interpreted sense) and use the fact the p-reflection is 
preserved by convergence in Hausdorff distance. 

□ 

Lemma A. 3. Let ilj for j G {1,2} and U be in Sr.R, o,i^d let a = a{r, R) e (0, 1) from the Holder 
estimates for harmonic functions in Lipschitz domains Lemma 3.13 Then the following estimates hold: 



\x[n,]-x[n2]\<r dff(r!i,r!2) 
||u[f]i]-M[n2]||ooS,K dH{ni,n2r 

\Jin^)~J{n2)\ <rM dH{^1.^2T 

|dist^(r2i,C/) -dist^(172,t/)| and |dist^(f/, r^i) - dist^(J7, l^a)! <r,fl 



(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 

(A.IO) 



Proof. We will start by showing (A.8), the only necessary assumption on the ^Ij is that one of their 
boundaries be rectifiable with -^^-i Hausdorff measure bounded in terms of r and R, 



< 



r,R 



dH{ni,n2) 



Then (A.IO I follows easily from (A.8), 

\(^t{U,VLi)-dM^{U,VL2)\ = 



d{x, dU)dx 



niiAu 

<R \niAn2\ 



n2Au 



d{x, dU)dx 



and because the dist is asymmetric we also check 



|dist {ni,U)- dist {n2,u)\ 



d{x, dfli)dx 



niiAu 

<R |f^lAf^2| + 

<R dHinun2) 
<R dHin,,n2) 



(OinOajACf 



d{x, dfl2)dx 

fl2AU 

\d{x, dfli) — d{x, dV,2)\dx 



Oinn2)AC/ 



dH{^i,^2)dx 



where in the third line we have used that if yi G dfli there exists 1/2 G dil2 with \yi — 2/2I < dni^i, ^2) 
and therefore d{x, 8^2) < d{x, dili) + dH{^i, ^2) and vice versa. Note that (A.9 1 is an easy consequence 
of ( [at] ) and ( |X6l ), 

Ui^i)- j{n2)\ < \\u[ni] - u[n2]\\^max{\ni\ , |f)2|} + l^^iArial, 



so we are left to prove the first three estimates. 

First recall the behavior of A under spatial dilations, 



X[an] = 



^ „-{N+2) 



x[n]. 



(All) 
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We will show that taking, 

one gets the following containments, 



ai7i C C —fli. 



(A.12) 
(A.13) 



This will imply by (A.ll) that 



and by rearranging, 



|A[r2i] - X[n2]\ < max{A[f7i],A[172]}(l 



a 



N+2 



Now to get an estimate of the form (A.6I, for t > note that 
to get 



|A[r!i 



1- + < {N + 2)t 

X[n2]\ < r-\N + 2)X[B,.]dH{dni,dn2). 



Let us prove (A. 131 for a chosen as in (A.12). If a; G ili \ let t e (0, 1) such that tx € dVL2. Then 
from the star-shaped property there is an exterior cone to Q.2 at te, 

te + C (x, 0t,) C \ r!2 where sinfl*, = j^, G (0, ^/2) 

so that 

t ' 

Then d{x, 8^1,2) > (i ~ l) ^^'^ moreover, 

1 



or by rearranging. 



- 1 r < dH{dQ.i,dVL2). 



a<t. 



Then by making the same argument for x £ Q2\^i we get 

ai7i C 172 C — f^i- 



Now we show (A.7|. Without loss of generality suppose that A[i7i] > A[il2]- Consider the difference 
of the droplet profiles on the two domains, 

w{x) — u[fli]{x) — u[il2]{x). 

First we estimate the size of w{x) on riiAi72. If a; G iliAil2 let h — d{x,d{^li U ^2)) then 
h <dH{^i,^2) and Bh{x) (lQiUn2 with dBh{x) n di^iU Q2) ^ ^■ 

By the standard construction of barriers for domains with the exterior cone property there exist < a < 1 
and C > depending on the uniform lower bound on the opening angle of the exterior cones such that 

\w{x)\<Ch- + ^h'. 
So on f2iAi72 and in particular on d{Qi D ^12) we have that 

\wix)\ S,fl max{ciff (f7i, ^2)", dff (Hi, 172)'} S,fl ^^(f^i, ^^2)", 
where the last inequality is due to the fact that (ijy(f7i, f72) < 2R. Finally on intersection x £ 17i n 172, 

Hx)\ <,R dH{n,,n2r + - H^2]\ <rM dHin,,n2r. 



Finally we show (A. 5 1. Let xq G 17i \ 172, then tx G M \ 172 for all i > 1. In particular let s > 1 such 



that xi — sxq G 917 1. We claim that: 



h d{xo, Q2) < d{xi,^2)- 



(A.14) 
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Let X e Bfi{xi), then consider the pomt s^^x, 

\s~^x — xo\ = s~^\x — xi\ < \x — xi\ = h. 
Therefore s~^x G Bh{xQ) C \ and so from the star-shapedness of SI2 the point x = ss~^x is as 



well, this proves the claimed inequality (A. 14). In particular: 



swp d{x,Q2) ~ sup d(x,dfl2)< sup d(x,dQ2) 
and by noting that fii and il2 play symmetric roles: 

dHinun2) <dH{dni,dn2). 



□ 



Define the metric space of boundaries of strongly star-shaped sets, 

dSr^R {dil : ft G Sr^n} 
with metric dH- This space embeds continuously into 

{/ € C°-\S^-') : r < / < i? and \\Df\\L^ < C{r,R)} 
with the L°° distance. As a direct consequence of this we get the following compactness for dSr.B- 

Lemma A. 4. (Compactness) The metric space {dSr.RT da) is compact: 

(i) Suppose that Tj e {dSr.Ri d-H) for some r,R> and all j e N. Then {Tj}j^fij has a subsequence 

that converges and any subsequential limit is also in dSr^R- 
(ii) Let I be a compact interval in K and Tj : I [dSr^R, dh) for j G N are an equicontinuous sequence 

of paths in {dSr^R, d^) ■ Then there is a subsequence of the that converges uniformly on I to 

a path r : / — > {dSr.R, dh). 

Proof. Arzela-Ascoli. The only interesting point is that dSr.R is closed in [JC^dn) where K. is the class 
of compact subsets of E^. It suffices to establish the estimate: 

dH{dni,dn2) <r,R dH{ni,n2) (A.15) 

where Vlj € Sr^R. Let x € dVli and y £ dCl2 such that: 

\x-y\^ d{x, dn2) = Sr^s). 

UxeR^\n2 then, 

dH{dQudn2) = \x-y\ = d{x, n2) < dH{fii,n2). 

Otherwise x e il2, then let t > such that (1 + t)x E dfl2- From the strong star-shapedness of fii we 
know that 

d{{l + t)x,ni)>t\x\f- > \{l + t)x~x\'- > Ix-yl"-. 

\x\ R R 

Rearranging we get that 

(aril, aria) ^\x-y\< -d{{i + t)x, n^) < -d^ (i^i, i^a). 

r r 

□ 
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Appendix B. Sup and Inf Convolutions 

For the comparison principle the convohitions in space only give spatial regularity of the free boundary 
at a point where a sup-convolved subsolution touches an inf-convolved supersolutions fr-om below. To fix 
this convolutions over space-time ellipsoids are used. These were introduced in ff 0' for viscosity solutions 



of the porous medium equation. Let u be a subsolution of (P-A| and v a supersolutionon M.^ x [0, -l-c»), 
let r > and c > 0. Define the sup-convolution of u 

u^'''^{x,t) = sup it(?/, s), (B.l) 

\{y,s)-{x,t)\<r-ct 



and the inf-convolution of ' 



inf viy,s). (B.2) 

\{v,s)-{x,t)\<r-ct 



For simplicity of presentation and since we want to emphasize the c — case separately we define 

u{x,t) = u'^'^{x,t) and v{x,t) — Vr,c{x,t). 

As for the convolutions in space only, c will boost the speed of the free boundary of v. 

Lemma B.l. (Space-time convolutions) 

(a) Letu as defined above the sup-convolution of u. 



(i) u is a subsolution of (P-A| in x [r, -|-cxi) with X{t) = sup|^_j|<^ A(s). 

(ii) The positive phase ^{u) has the interior ball property. More specifically, let Pq — (xo^to) G 
dfl{u) then the positivity set of u has an exterior ball of radius r centered at Pq, 

SriPo) {(y.s) e X [0,+oo) : \{y,s) - (a;o,to)| < r} C n{uf 

and there exists Pi G dBr{Po) H 90 (u). Then Br (Pi) touches 17+ (w) from the inside at Pq. 
(b) Letv as defined above the inf-convolution of v. 

(i) v is a supersolution of ( |P-A ) in 'U.^ x [r, r/c) with \{t) — inf|s_t|<,._ct ^{s). 



(ii) The positive phase ^{v) has an exterior ellipsoid at every point of its free boundary. More 
specifically, let Pq G dQ-i-{v) then Br-ctoi-Po) touches the free boundary of v from the inside i.e. 

Br-cto{Po) ^ ^{v) and there exists Pi G dBr-cto{Po) ^ dfl{v). 
The ellipsoid 

Er,s{Pi) = {(y-s) : Ky.s) - <r~cs} 
touches from the outside at Pq. 
(Hi) Let Hi be the hyperplane orthogonal to Pi ^ Pq passing through Pi, i.e. the tangent hyperplane 
to Br-ctoiPo) o,t Pi o-iT-d- {i^,m) G M" x be the normal vector to Hi rescaled so that = 1. 
We call m the advance speed of the free boundary of v at Pi. This terminology is justified when 
V G C'^'^iW^ X M+). Let Hq be the tangent hyperplane to Er^s{Pi) cit the point Pq and (u' ,m') 
it's normal vector rescaled as before. Then m' = m-\-5, the inf-convolution boosts the advance 
speed of the free boundary by 5. 

Appendix C. Higher Regularity for the Equilibrium Problem 



In this appendix we describe the method in 13 and 22 to show that if u is a solution of the 
equilibrium problem 

-Aw(a;) = A in f7 = VL{u), 

dEQ) 



\Du\ = 1 on r = T{u). 



with n G 5o then F is C^-". The first step is Theorem [2^ from |i3] which says that F is C^'". As 
a consequence of this u G C^ °'{Vt). Then one uses a transformation of the problem first used for this 
purpose by [22] that has the effect of flattening the free boundary while making the PDE and the 
boundary condition nonlinear. Below we describe the transformation in detail. 
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We may assume that £ F and 1 = |I?u|(0) = 9„u(0). Then near we make the change of variables: 
Vj^Xj for l<j<n-l, Vn^u, w = x„, 
Then we can write the y derivatives of v in terms of the x derivatives of u so that v{y) solves the problem: 

\lT,=id,{F,{Dv)) = \ for y„>0 



{dnvf - ELi {djvf = 1 on 2/„ = 



where 




Pn 



for 1 < i < ~ 1, 



(C.l) 



(C.2) 



{Pn? 

It is straightforward to check that this operator is uniformly elliptic when dnV is bounded from below. 



This fact as well as the nondegeneracy of dnV from the boundary condition enables the argument in 22 
to go through for our problem and yield the desired regularity result in Corollary |2.3| 



Let us show that (C.l) is the correct equation. Based on the definitions of y.j and w, we have 

u(yi,...,?/„_i,w(yi,.. .,?/„)) = y„. (C.3) 



Taking the n-th partials in (C.3) we get 

1 



dnnV 



and dnj u 



dnnVdjV 



dnjV 



where any derivative of u is assumed to be evaluated at (yi, yn-i, v{yi, ...,?;„)) and any derivative of 
V at (j/i, 2/„). Similarly taking derivatives in (C.3| for 1 < j < n — 1, we have djU + dnudjV = and 



thus 



'djjU ~ 2dnjUdjV + dnnU{djvY + dnUdjjV 



[dnVY dnV 



Then we substitute into (EQ) to get 



n-l 

A = - ^ djjU - 

ri-l 



n-l 

E 



n-l 



dr, 



{dnVf 



J = l 
n-l 



{dnVf 



djV ' 

dnV, 



- dn 



2 {dnVY 



The calculation of the boundary condition is easier 



1 = \Du\' = 



n-l 

E 



d„v 
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